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ABSTRACT 

The Panchromatic Hubble Andromeda Treasury (PHAT) is an on-going Hubble Space Telescope 
(HST) Multicycle Treasury program to image ~l/3 of M31's star forming disk in six filters, spanning 
from the ultraviolet (UV) to the near-infrared (NIR). We use the Wide Field Camera 3 (WFC3) 
and Advanced Camera for Surveys (ACS) to resolve the galaxy into millions of individual stars with 
projected radii from 0-20 kpc. The full survey will cover a contiguous 0.5 square degree area in 828 
orbits. Imaging is being obtained in the F275W and F336W filters on the WFC3/UVIS camera, 
F475W and F814W on ACS/WFC, and F110W and F160W on WFC3/IR. The resulting wavelength 
coverage gives excellent constraints on stellar temperature, bolometric luminosity, and extinction for 
most spectral types. The data produce photometry with a signal-to- noise ratio of 4 at rap275W = 25.1, 
^F336W = 24.9, ra F 475W = 27.9, ra F 8i4W = 27.1, ra F iiow = 25.5, and ra F i60W = 24.6 for single 
pointings in the uncrowded outer disk; in the inner disk, however, the optical and NIR data are 
crowding limited, and the deepest reliable magnitudes are up to 5 magnitudes brighter. Observations 
are carried out in two orbits per pointing, split between WFC3/UVIS and WFC3/IR cameras in 
primary mode, with ACS/WFC run in parallel. All pointings are dithered to produce Nyquist- 
sampled images in F475W, F814W, and F160W. We describe the observing strategy, photometry, 
astrometry, and data products available for the survey, along with extensive testing of photometric 
stability, crowding errors, spatially-dependent photometric biases, and telescope pointing control. We 
also report on initial fits to the structure of M31's disk, derived from the density of red giant branch 
stars, in a way that is independent of assumed mass-to-light ratios and is robust to variations in dust 
extinction. These fits also show that the 10 kpc ring is not just a region of enhanced recent star 
formation, but is instead a dynamical structure containing a significant overdensity of stars with ages 
>lGyr. 

Subject headings: galaxies: stellar content — stars: general — stars: imaging — galaxies: individual 
(M31) 
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1. INTRODUCTION 

Our quest to understand the Universe relies on detailed 
knowledge of physical processes that can only be cali- 
brated nearby. It is impossible to interpret observations 
across cosmic time without an underlying understand- 
ing of stellar evolution, star formation, the initial mass 
function, the extinction law, and the distance scale, all 
of which require detailed studies of individual stars and 
the interstellar medium (ISM) on sub-kiloparsec scales. 

When the needed studies of stars and gas are carried 
out in the Milky Way, they frequently face complica- 
tions from line-of-sight reddening, uncertain distances, 
and background/foreground confusion. As such, it is 
sometimes easier to constrain physical processes in ex- 
ternal galaxies, which are free of the projection effects 
that can plague Milky Way studies. Not only are obser- 
vations of external galaxies more straightforward to in- 
terpret, but they can also be placed in the larger context 
of the surrounding environment (i.e., the ISM, metal- 
licity, and star formation rate (SFR)). Galaxies in the 
Local Group therefore offer an excellent compromise be- 
tween being close enough to resolve relatively faint stars, 
while being distant enough to unveil the complex pro- 
cesses that govern star and galaxy evolution in their full 
galactic context. 

Unfortunately, even the nearest massive galaxies have 
sufficiently high stellar surface densities that severe 
crowding compromises the detection of fainter, more age- 
sensitive stellar populations, allowing only the brightest 
stars to be studied with typical ground-based angular res- 
oluti on in high-surface br ightness regions of galaxy disks 
(e.g., [Massey et aL]|2QQ6[ ). However, with the high angu- 
lar resolution available from HST, we have the potential 
to resolve millions to billions of stars within the Local 
Group, grouped into galaxies with a common distance 
and foreground extinction. These stars, along with their 
ancestors and descendants (e.g., molecular clouds, Hll re- 
gions, variable stars, X-ray binaries, supernova remnants, 
etc.), provide transformative tools for strengthening the 
foundation on which knowledge of the distant Universe 
is based. 

Within the Local Group, the Andromeda Galaxy 
(M31) offers the best proxy for the properties of more 
distant galaxies. It is massive (sampling above the char- 
acteristic stellar mass (3 — 5xlO lo M ) over which rapid 
systematic changes in galaxies' stellar pop ulations and 
structure occur (e.g., |Kauffmann et al.|2003| )), hosts spi- 
ral structure, and contains the nearest example of a tra- 
ditional spheroidal component (outside the MW). M31 is 
also representative of the environments in which typical 
stars are found today. More than half of all stars are 
currently found in t he disk and bulges o f disk-dominated 
galaxies like M31 (Driver et al. 2007), and more than 
3/4 of all stars in trie Universe have metall icities within 
a factor of two of solar ( jGallazzi et al.|2 008 ) , comparable 
to the typical metallicities of stars in M31. 

In addition to its dominant solar-metallicity popula- 
tion of young stars, M31 also contains significant popu- 
lations of older super-solar metallicity stars in its bulge, 
and sub-solar metalli c ity stars in its outskirts (e.g., 
Sirchenko et al.||1998l |Lauer et al. | [20121 JBrown et al. 




making M31 a superb laboratory for constraining stellar 
evolution models across at least an order of magnitudue 
in iron abundance. Moreover, because of its high mass, 
M31 contains >90% of the stars in the Local Group (out- 
side the MW), making it ideal for generating samples of 
sufficient size that Poisson statistics are negligible and 
even rare phenomena are well represented. Finally, the 
stars in M31 are bright enough to be accessible spectro- 
scopically, allowing one to augment imaging observations 
with spectroscopy, providing measurements of the kine- 
matics, metallicities, spectral types, and physical param- 
eters of star clusters and massive main sequence, asymp- 
totic giant branch (AGB), and red giant branch (RGB) 
stars. 

We have therefore undertaken a new imaging survey 
of M31's bulge and disk using HST. The Panchromatic 
Hubble Andromeda Treasury (PHAT) survey is being 
carried out as a "Mult i- Cycle Treasury" program to im- 
age a large contiguous area in M31, building upon the 
existing ground-based studie s that probe M31's most lu- 
minous i 
|etH 

survey uses HSTs new instrumentation to provide spec 
tral coverage from the UV through the NIR, with which 
one can effectively measure the bolometric luminosity, 
spectral energy distribution, and morphology of most as- 
trophysical objects. Broad-band coverage allows one to 
constrain the mass, metallicity, and ages of stars, even in 
the presence of extinction. Other objects such as back- 
ground active galactic nuclei (AGN), planetary nebulae, 
X-ray binaries, and supernova remnants, whose bolomet- 
ric luminosity peaks outside the accessible wavelength 
range, will still have distinctive spectral energy distribu- 
tions within the HST filter coverage. Thus, the broad 
wavelength coverage of the survey enables full charac- 
terization of stars, their evolved descendants, and useful 
background sources. When complete, the survey should 
contain data on more than 100 million stars, comparable 
to the number found in the Sloan Digital Sky Survey. 

The legacy and scientific value of our survey is rich and 
diverse. The survey data can be used to: (1) provide the 
tightest constraints to date on the slope of the stellar 
initial mass function (IMF) above > 5 M as a func- 
tion of environment and metallicity; (2) provide a rich 
collection of clusters spanning wide ranges of age and 
metallicity, for calibrating models of cluster and stellar 
evolution (e.g., I Johnson et al.|2012| ); (3) characterize the 
history of star formation as a function of radius and az- 
imuth, revealing the spiral dynamics, the growth of the 
galaxy disk and spheroid in action, and the role of tidal 
interactions and stellar accretion; (4) create spatially 
resolved UV-through-NIR spectral energy distributions 
of thousands of previously cataloged X-ray binaries, SN 
remnants, Cepheids, planetary nebulae, and Wolf-Rayet 
stars, allowing full characterization of these s ources; (5) 



2003T Worthey et al.| 120051 IKalirai et al, 



2006a| |Chap 



man et al.||2006| Fan et al.||2008[ and references therein) , 



constrain r are phases of stellar evolution (e.g., Rosenfield 
|et al.|2012| ); and (6) provide rich probes of the gas phase 
and its interaction with star formation, by using sub- 
arcsecond extinction mapping as a probe of the molecular 
gas, and by comparing parsec-scale gas structures to the 
recent history of star formation and stellar mass loss. In 
the future, the survey will provide the fundamental base- 
line for characterizing sources in transient surveys (e.g., 
PanSTARRS, LSST) both by using direct identification 
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of counterparts, and by associating transients with the 
properties of the surrounding stellar population, when 
the counterpart is undetected even at HST's depth and 
resolution. The survey will also catalog hundreds of UV- 
luminous background AGN, which can be used as ab- 
sorption line probes of the detailed physics of the ISM in 
future UV spectroscopic observations, and as a reference 
frame for proper motion studies. 

In this paper after giving a brief history of stellar pop- 
ulation studies in M31 (Section|2]), we present the survey 
design (Section [3| and data reduction (Section [4|, along 
with a thorough characterization of the data quality. In 
Section [5j we present initial color magnitude diagrams 
covering several kiloparsec-scale contiguous regions, lo- 
cated in the bulge, two major star forming rings, and 
the outer disk. We also present a preliminary analysis 
of the structure of M31's disk, based on counts of RGB 
stars in the NIR. 

2. THE HISTORY OF M31 RESOLVED STELLAR 
POPULATION STUDIES 

The recognition of the prese nce of di f ferent stellar pop- 
ulations in M31 began when iHubble (fl929|) noted that 



([192! 

esolv 



the outer parts of the disk could be resolved into stars, 
while the central region showed only unresolved light, de- 
spite clear spectroscopic evidence that it was made up of 
stars. When Baade carried out his 20-yea r study of M31 
(see |Baade fc Payne-Gaposchkin||1963| ), he used red- 
sensitive plates to show that the central regions could 
indeed be resolved into stars, and had properties similar 
to those of Galactic globular clusters. Based on these 
M31 studies, Baaded developed the influential concept 
of Populations I and II, which formed the foundation 
of subsequent population studies. The central areas of 
M31 were made up of red, low-luminosity stars (Popula- 
tion II) and the main disk contained luminous blue stars 
(Population I). At the time of Baade 's work, theor etical 
stellar evolution calculations ( Schwarzschild 1965) were 
beginning to provide models that nicely tit the kinds of 
patterns found observationally by Baade, allowing the 
empirical facts of different stellar properties to be turned 
into quantitative population histories. Baade concluded 
that the central bulge of M31 was made up of very old 
stars, while the bright blue arms were young. He mapped 
out the structure of these arms and showed that the most 
luminous areas were fragments of arms located at r adial 



distances of between 8 and 12kpc (see also work by Arp 
1964). Subsequent work by m any astronomers (see 



van 



den bergh 1991 Hodge 1992 for references) confirmed 
this pattern, but showed that the dichotomy of just two 
populations was too simp l e. Po pulation studies of the 
disk of M31 by Will iams] ( |2QQ3[ ), for example, showed 
that the populations vary across the disk, implying dif- 
ferent histories of star formation within the broader clas- 
sification of Population I. The structures seen between 8 
and 12kpc are now thought to be a ring of star forma - 
tion at - 10 kpc flHabing et al.|1984| |Gordon et al.|2006[ ) , 
based in large part on observations of M3Ts interstellar 
medium. 

Global population studies were sparse in the later parts 
of the 20th century, in spite of major wide- area ground- 
based CCD surveys o f M31's bright stellar content (e.g., 
Magnier et al. 1992). Most papers instea d looked at 
smaller units, e.g., globular clusters (e.g., Rich et 



1996 |, OB associations (e.g., jVlassey et al.|1986| He 
et al.J1994l|Hunter et al.| lg96p , and individual HST 



Haiman 



3.g., [Rich fc Mighell| |l995 ) . 

iry, large new catalogs of M 



tury 

the optical dMassey et al 
et al.j2006| [Mould et al.|2UU8p 



and individual HST fields 
After the turn of the cen- 
31 's stars were pr oduced in 
20061) and infrared 



(Skrutskie 



but much of the more re- 



cent work on stellar populations c oncentrated on the halo 
and on extended disk stars (e.g., [C uillandre et al. 2001; 
Durrell et al.|200l||Ferguson fc Johnson|2001||!Sa raj edini 
fc Van Duyne||2001l|Rich et al.||2004UBrown etaIT2 006 ; 
Kalirai et al.||2006a|bl |Brown et al. | |20071 120081 |Richard-| 
son et al.||20081 |Brown et al.||2009b| |Bernard et al.||2011[ 
and many others), and the popula t ion of the bulge (e.g 



Saglia et al.l[2010l [Davidge 2001 
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2005 

Sarajedini fc Jablonka 2005; Stephens et al. 2003; Olsen 
jet al ||2006| |Rosenfiel d et al.||2012p , which holds special 
interest due to the presence of a significant super-solar 
meta llicity stellar population. Analyses of stellar clusters 



, Krienke fc H odge 2007, 2008; Barm by et al.||2009 
Hodge et al. 2009; Perina et al. 2Quj|) a nd various HS 



point ings have also continued in earnest (Bellazzini et al. 



2003J). 

These and other studies have confirmed a basic pic- 
ture where M31 hosts a clear disk and bulge. How- 
ever, there is evidence for more complex structures in 
the inner region (including a bar and a boxy pe anut- 
shaped bu l ge, in addition to M31's classical bulge; [ Lind- 
blad 19561 |Stark|[l977l IStark fc Binney||19941 |Athanas- 



|soula fc Beaton|2006||Beaton et al.|2007p , the disk (which 
shows a change in p osition angle at ^18 kpc, most likely 
Walterbos & Kennicutt 1988), and 

1 ? : r-i I m i II i ' i T 



due to a warp; e.g 
the outer disk and halo 



most recently, Tempel et al 



2010 McConnachie et al. 20091). The disk contains am- 
ple evidence tor recent star formation, confined largely 
to major spiral arms or the 10 kpc ring. Near the ring, 
M31's current metallicity appears to be comparable to 
or higher than that of the Milky Way, with ambiguous 
evidence for a gradient in O/H; however, existing metal- 
licity data is surprisingly sparse, spanning a limited range 
in radii with larg e 0.5 dex scatter i n the metallicity at 
any given radius dRubin et al.||1972| |Dennefeld & Kunth 



1981 Bl air et al.| 1982; Zarits ky et al.|1994||Galarza e~ 



1999 



2001 



IVenn et al.||2000| ISmartt et al ||2001| |Han etTal 
Trun dkTet al.|2QQ2 ). Analyses of older stellar pop- 



ulations (RGB and PNej suggests that the dis k of M31 
hosts a wide range of stellar metallic i ties (e.g., Wortheyj 
eTaL||2005l |Jacoby fc Ciardullo|[l999l |Richer et al"1f l999 ] 
and references therein). The metallicity of older stars 
in M31's bulge likewise spans a wide range, but seems 
to reach super-solar metallicities in the v ery center (e.g., 



Sirchenko et al.p998| |Lauer et al.||2012l and references 
therein) . 

Somewhat surprisingly, the distance to M31 remained 
largely uncontroversial in the last two decades, unlike for 
the other galaxies in t he Local Group, be it the Large 



Magellanic Cloud (e.g., Macri et al. 2006; [van Leeuwen 
l et al.|2007l) or M33 (e.g., |Bonanos et al ||2U06[ |Scowcrott = 
|et al.|20(J!?fl . For example, th e Ceph eid dist ance mo dulus 
to M3T ot \l = 24.41 ± .08 ( |Freedman fc Madore|1990| 
see also |Riess et al.|2012[ using data from the PHAT sur- 
vey) agrees well with the red clump distance modulus of 
/i = 24.47±0.06 ( |Stanek fc Garnavich|1 998[), a value nu- 

ased 



merically ide ntical 



to the TR GB- based distance of fiQ = 



al] 24.47 ± 0.07 ( |McConnachie et al.|2Q05| and very close to 
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the /io = 24.44 ± 0.12 measure d via a direct met hod us- 
ing a detached eclipsing binary ( Ribas et al.|2QQ5 ). Given 
such too-good-to-be-true agreement one would normally 
suspect a "bandwagon effect", but all these papers use 
different methods and different zero-point calibrations to 
derive their distances. We will therefore adopt an M31 
distance modulus of /im 31,0 — 24.45 ±0.05 (physical dis- 
tance of dj\43i = 776 ± 18 kpc) in this paper. 

3. SURVEY DESIGN 

The optimal survey design for efficiently imaging a 
large portion of M31 requires finding solutions to several 
problems. These include determining the best portion of 
the galaxy to image, the best filters to use, and the best 
design of the exposure sequences. These in turn require 
identifying the best way to efficiently tile large areas with 
two cameras operating in parallel, while maximizing the 
photometric depth and image resolution, given a limited 
total exposure time for any location. In this section we 
discuss the rationale for the choices we made to address 
all of these issues in a way that we hope has optimized 
the broad scientific u tility of the program. 

We start in Section I37T1 with the rationale for the loca- 
tion and geometry of the PHAT survey area. Section |X2 
gives the scientific rationale for 6-filter coverage from the 
UV through th e N IR, and describes the specific filter 
choices. Section pO] explains the exposure strategy for fit- 
ting observations in six filters into two orbits. Section 3.4 
shows how the observations are dithered to optim ize s pa- 
tial resolution and repair detector defects. Sect ion [33] de- 
scribes how the 2-orbit visits are packed into an efficient 
3x6 mosaic of 18 pointings, building the 23 " bricks" that 
tile the survey area. We close in Section 3.6 with a brief 
discussion of a coordinated spectroscopic campaign. 

3.1. Areal Coverage 

We evaluated several possible schemes for mapping 
M31. We initially considered producing a complete map 
of M31's star forming disk, using only one orbit per point- 
ing and a reduced number of filters. However, meeting 
the broadest set of scientific goals required more com- 
plete filter coverage, and thus longer exposures at each 
position. We therefore reduced our areal coverage to a 
generous quadrant of M31. Since spirals like M31 have 
fairly regular structure, we judged that a quadrant would 
be sufficient to characterize the galaxy while still cover- 
ing enough area to allow statistically significant samples 
of previously cataloged objects, of star forming regions 
across the widest possible range of star formation in- 
tensity and metallicity, and of azimuthal variations in 
the star formation rate. We specifically chose the north- 
east quadrant, which has the lowest internal extinction, 
the largest number of regions with unobscured, high- 
intensity star formation, and the least contamination 
from M32. However, the tiling is generous enough that a 
substantial fraction of the northwest quadrant is covered 
by the HST imaging as well. 

The resulting survey area is a long, roughly rectangular 
region with a slight bend in the middle. The long axis 
of the survey begins at the center of the galaxy, covering 
much of the bulge, and extends to the northeast following 
the major axis of the disk out to the last obvious regions 
of star formation v isible in GALEX imaging (see maps in 
Thilker et al.|2005[ ) . The short axis spans the minor axis 



of the galaxy, out to a comparable projected radius as the 
long axis, once M31's inclination is taken into account. 

In richer detail, the survey area is contiguous and is 
divided into 23 sub-areas, called "bricks," which solve a 
number of problems related to the efficient design and 
operation of a large multi-year surv ey; we discuss the 
design of the bricks in detail in section 3.5 In the present 
context they can be considered to be rectangular tiles 
that "pave" the complete survey area. The bricks are 
arranged into two strips that respectively comprise the 
northwest and southeast halves of the survey area. Their 
detailed orientation and positioning gives the edges of the 
survey area a somewhat serrated appearance. Maps of 
the brick positions are shown in Figure [T] superimposed 
on various mult i- wavelength images. 

The northern strip has 12 bricks extending along the 
major axis from the bulge to the outer disk. The bricks 
in the southern strip are aligned to complete coverage 
of the quadrant. The naming scheme is such that Brick 

1 covers the bulge, with odd numbered bricks extending 
along the major axis. To increase the windows in which 
observations can be scheduled, Bricks 1-11 are observed 
with a slightly different orientation than Bricks 12-23. 
In practice, we use the brick designations as a shorthand 
for referring to different regions of interest within the 
survey. 

The survey area samples the diversity of environments 
within M31. The southwest end of the survey encloses 
much of the bulge, including the nucleus and the dense 
stellar population that surrounds it. The extent of the 
survey along both the major and minor axes of M31 is 
sufficient to observe the transition of the bulge into the 
main disk of M31. Star formation is probed at many 
locations within the disk, the most notable being the 
10 kpc ring of strong star formation, which is tracked 
throughout the southern edge of the survey. The width 
of the survey area completely tracks the roll-off of this 
arm into the background disk, and on the southern side 
of the galaxy, the terminus of the star forming disk itself, 
into the surrounding halo. Weaker arms of star formation 
just outside the bulge are sampled, as are weaker spiral 
arms well outside of the main ring of star formation. 
Between the zones of strong star-formation, the survey 
samples the smooth background population of the disk 
over the entire major axis of M31. The northeast end of 
the survey also captures the termination of the disk and 
the transition into the halo. To a fair approximation, the 
bricks that tile the survey area typically sample one or 
two of these various features. Table [T] presents a brief 
description of the populations present in each brick. 

Because of the large number of orbits required by this 
program, not all bricks can be observed in a single Cycle. 
We have therefore prioritized the bricks to maximize the 
possible science output in early years. Our Year 1 pri- 
ority was to complete bricks that sample the full radial 
extent of the galaxy along the major axis, focusing on the 
bulge and major star forming rings and/or spiral arms 
(Bricks 1, 9, 15, 21, followed by Bricks 17 and 23). Year 

2 priorities were to sample the high intensity star form- 
ing ring (Bricks 2, 8, 12, 14, and 22), and to increase 
the sampling of star forming regions on the major axis 
(Brick 5). Year 3 priorities are to complete most obser- 
vations of the major star forming ring (Bricks 4, 6, 16, 
and 18), and to start building larger contiguous regions 
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in the inner and outer galaxy (Bricks 3, 19). Year 4 will 
be devoted to completing all remaining areas (Bricks 7, 
10, 11, 13, and 20). The status of observations (as of Fall 
2011) are given in Table [T] 

3.2. Filter Choices 

The choice of UV-through-IR filter coverage was driven 
by a number of goals. The need for two optical filters 
was obvious, given that optical HST imaging data has 
proved to be the most efficient route to deriving the 
deepest possible color-magnitude diagrams (CMDs) for 
the largest number of stars. Supplementing the optical 
data with two additional NIR filters allows one to extend 
stellar population studies to dusty regions, and to bet- 
ter constrain the bolometric fluxes of intrinsically cool 
stars in important evolutionary phases (AGB stars, Car- 
bon stars, red core Helium-burning stars). Adding an 
additional two UV filters opens up science that can be 
done with hot stars, and in particular permits simultane- 
ous constraints of effective temperature and extinction, 
when combined with measurements in optical filters ( e.g. , 



Bianchi|2007| |Romaniello et al.||2002| |Zarits ky 1999 ) . 
The choice oi specific filters, shown in Figure [2j \ 



was 

made as follows: 

In the UV, we adopted F336W as the filter giving the 
highest throughput. This filter is also immediately blue- 
ward of the Balmer break, giving the best constraint on 
its amplitude. For the bluer UV filter, we wanted to 
push to the shortest possible wavelengths to give the best 
constraints on stellar temperature for the hottest stars, 
and to provide the largest baseline for constraining the 
extinction law. However, the bluest wide WFC3/UVIS 
filters have relatively low efficiency, and moreover fall 
within the 2 175 A dust extinction feature, which is known 
to vary dramatically with environment. We therefore 
adopted F275W as the bluest, high throughput UV filter 
that would not be significantly affected by variations in 
dust composition. 

In the optical, we adopted the F814W filter, which 
has consistently offered the highest throughput for stel- 
lar populations studies. The choice of the bluer optical 
filter was more problematic, as several valid choices exist. 
The F606W filter offers very high throughput and has 
been widely used for previous stellar population stud- 
ies. However, it is quite red, providing a limited color 
baseline in the optical, and weaker constraints on the 
amplitude of the Balmer break (when combined with 
F336W). The F555W filter provides a wider color base- 
line than F606W, and has also been used for a large 
number of HST stellar population studies. However, it is 
significantly narrower, and thus has much lower through- 
put. We therefore rejected these two filters and adopted 
the F475W filter (approximately SDSS-#). This filter 
is as broad as F606W, but is bluer than both F606W 
and F555W, providing much better color separation from 
F814W. When combined with F336W, it provides good 
constraints on the Balmer break. The only cost is some 
loss in depth for intrinsically red stars (i.e., on the red gi- 
ant branch). However, since much of the disk is crowding 
limited, this limitation is not severe at most pointings. 

For the NIR, we use the F110W and F160W filters. 
These two filters provide the highest throughput with the 
WFC3/IR camera, and have been used successfully in our 
previous SNAP survey of stellar populations in nearby 



galaxies ( |Dalcanton et al. 2012 ) . The only drawback with 
this combination is the partial overlap of the F110W and 
F814W bandpasses. The only other feasible substitutes 
for the F110W filter would have been the F140W filter, 
which overlaps the F160W filter, and the F125W filter, 
which is much narrower than F110W, and which has less 
color separation and temperature sensitivity when paired 
with F160W. 

Almost all of the PHAT filters have been used in cal- 
ibration observations of nearby stellar clusters. The 
WFC3 Galactic Bulge T reasury Program (GO-11664; 
|Brown et al.||2009a| |2010|) have taken calibration obser- 
vations in F814W, ihlOW, and F160W. F336W obser- 
vations of the same clusters have been carried out in 
GO- 11729 (PI: Holtzman). Our survey is executing ob- 
servations of the same clusters in combinations of F275W 
and F475W filters, to complete the calibrations of our fil- 
ter set. Additional calibration in F336W, F475W, and 
F814W will be provided by GO-12257 (PI: Girardi) for 
intermediate-age Magellanic Cloud clusters. 

The adopted filter set should allow us to make strong 
constraints on the effective temperature and extinction 
of the stars in our sample. With the UV+optical fil- 
ters, we expect to be able to separate T e ff and E(B— V) 
with little degeneracy for both hot stars (12,000 <T e /-/- < 
40,000if; IMassey et aLl[l995l jRomaniello et al.|[2 Q02; 



| Bianchi|[2007p and cooler stars (5500-6500/f ; |Zaritsky| 
1999). Inspection of reddening- free diagrams suggests 
that the optical+NIR combination will allow us to ex- 
tend E(B— V) constraints to cooler stars (T e ff < 5000K) 
as well. 

3.3. Exposure Sequences 

The primary aims of our exposure plan are (1) imag- 
ing two filters per camera; (2) achieving Nyquist-sampled 
images through dithering where possible; and (3) avoid- 
ing saturation of bright sources. As we describe below, 
these goals are balanced against constraints on the num- 
ber of images that can be downloaded when running 
WFC3 and ACS in parallel, and on limitations on the 
duration of an orbit. Because of the strains that this 
program puts on HST's observing schedule, the expo- 
sure sequences must fit within the shortest possible orbit 
duration ("schedlOO", i.e., an assumed duration that is 
schedulable for 100% of the orbits). This constraint max- 
imizes the schedulability . par ticularly in the summer ob- 
serving season (Section |3.5[ ) which has a more restric- 
tive scheduling window. In the winter observing season, 
we can relax the orbit length constraints to "sched60" 
(a duration that fits within 60% of the orbits), giving 
slightly longer exposures. 

Scheduling observations in six filters requires 
two orbits, with the first orbit devoted to 
WFC3 /UVIS+ ACS /WFC and the second to 
WFC3/IR+ACS/WFC; this ordering minimizes persis- 
tence in the WFC3/IR channel (ISR WFC3 2008-33; 
McCullough & Deustua 2010), by allowing more time 
for the persistent charge to decay. We run WFC3 in 
primary, and ACS in parallel^ 

21 Because the corrections for differential velocity aberration 
(ISR OSG-CAL-97-06; Cox 1997) by the pointing control software 
are optimized for the primary observations, the ACS parallels will 
not automatically be corrected for this milli-arcsecond level effect. 
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The strongest constraints on the observing sequence 
within each orbit come from the limited buffer space 
available on-board the HST instruments. The ACS 
buffer can only hold one full- frame ACS/WFC image, 
and the WFC3 buffer can hold either two full-frame 
WFC3/UVIS images, or three WFC3/IR images (if the 
number of no n- destructive reads in the latter is reduced 
to ~10 from the nominal 15 samples) before needing to 
be dumped. However, the time to clump the buffer is 
substantial (~340 seconds), which can lead to significant 
inefficiencies if the observing sequence does not allow the 
buffers to dump in parallel with the observations. This 
issue is particularly severe when running both imaging 
cameras simultaneously, because of the high data rates. 

After extensive experimentation, we were able to find 
exposure sequences that allow four WFC3/UVIS expo- 
sures and four ACS/WFC exposures in the first orbit, 
and five WFC3/IR exposures and five ACS/WFC expo- 
sures in the second orbit, with minimal latencies due to 
buffer dumps. The resulting observing sequence also fits 
in dithers between every pair of exposures, with the ex- 
ception of the first two in the WFC3/UVIS+ ACS/WFC 
orbit (see Figure [3]). This high observing efficiency comes 
at the expense of unequal exposure times for observa- 
tions in a single filter, and fewer non-destructive reads in 
the WFC3/IR observations (to allow the buffer to hold 
more WFC3/IR images between dumps). The sequence 
is summarized in Table [2j 

During the first orbit, the four WFC3/UVIS exposures 
follow the sequence: F336W (550 sec; 550 sec), F275W 
(350 sec; 350 sec), F336W (700 sec; 800 sec), F275W 
(575 sec; 660 sec), where the two numbers in parenthe- 
ses indicate the exposure times in the summer schedlOO 
and winter sched60 observing seasons, respectively. We 
obtained two exposures in each filter to enable CR re- 
jection and cover the chip gap; three exposures per filter 
would have been preferable, but constraints on buffer 
dumps limited the observations to no more than four 
WFC3/UVIS exposure per orbit. Note that the major- 
ity of the chip gap is only imaged in one exposure per 
filter, making CR rejection more challenging in this re- 
gion. 

The four ACS/WFC exposures in the first orbit 
are all in F814W, with exposure times of (15 sec; 
15 sec), (350 sec; 350 sec), (700 sec; 800 sec), and (455 sec; 
550 sec). The short 15 sec F814W exposure is included to 
allow photometry for stars brighter than F814W^ 17. 5p] 

During the second orbit, the five WFC3/IR expo- 
sures follow the sequence F160W (NSAMP=9, STEP200, 
399 sec; NSAMP=9, STEP200, 399 sec), F110W (NSAMP= 
13, STEP100, 699 sec; NSAMP=11, STEP200, 799 sec), 
F160W (NSAMP=9, STEP200, 399 sec; NSAMP=9, STEP200, 
399 sec), F160W (NSAMP=9, STEP200, 399 sec; NSAMP=9, 
STEP200, 399 sec), F160W (NSAMP= 9, STEP200, 399 sec; 

The affine corrections used for astrometry (Section |4.6| > should pro- 
vide adequate corrections, however. 

22 There is no equivalent short "guard" exposure in the UV, 
since the number of saturated stars was expected to be n egligible 
( based on existing ground-based {/-band data from the |Massey| 
et al. ( 2006 ) Local Group Survey), and the penalty in exposure time 
would be large for the unsaturated stars (given the small number of 
possible WFC3/UVIS exposures). In practice, it may be possible 
to pull out reasonable photometry from even the saturated stars 
(see ISR WFC3 2010-10 b y Gilliand et al. 2010 for WFC3/UVIS 
and |Anderson et al.|2008| for ACS/WFC). 



NSAMP= 11, STEP100, 499 sec), where the numbers in 
parentheses give the number of samples, the MULTIACCUM 
exposure sequence, and resulting exposure time for the 
summer schedlOO and winter sched60 observing sea- 
sons, respectively. The STEP sequence was adopted to 
provide optimal sampling for the wide range of fluxes ex- 
pected for observations of stellar populations. Since the 
MULTIACCUM sequence allows for CR rejection in an ac- 
cumulated single image, multiple images are not needed. 
We therefore chose to use multiple pointings in F160W 
only, to allow the construction of a Nyquist-sampled im- 
age in this filter, at the expense of not being able to reject 
chip defects in F110W (the WFC3/IR "blobs" and the 
"death star"; ISR WFC3 2010-06, Pirzkal 2010). This 
trade-off was preferred, given the lower resolution of the 
WFC3/IR channel compared to the other cameras, and 
the partial overlap of the F110W and F814W bandpasses. 

The five ACS/WFC exposures in the second orbit are 
all in F475W, with exposure times of (10 sec; 10 sec), 
(600 sec; 700 sec), (370 sec; 360 sec), (370 sec; 360 sec), 
and (370 sec; 470 sec). As for F814W, we include a short 
"guard" exposure to allow photometry of stars brighter 
than F475W- 18.5. 

The resulting total exposure times for each two- 
orbit visit are (925 sec; 1010 sec) in F275W, (1250 sec; 
1350 sec) in F336W, (1720 sec; 1900 sec) in F475W, 
(1520 sec; 1715 sec) in F814W, (699 sec; 799 sec) in 
F110W, and (1596 sec; 1696 sec) in F160W, where the 
first and second numbers indicate the exposure times in 
the summer and winter observing seasons. The effective 
exposure times for the ACS observations will be at least a 
factor of two larger due to overlaps (see Figure pi), giving 
nearly two full orbits of exposure time in botn F475W 
and F814W outside the chip gap. 

3.4. Dithering Strategy 

The exposure sequences were interleaved with small 
angle maneuvers to produce dithered images. Large 
dithers provide for the repair of detector defects (hot 
or bad pixels, missing columns, and so on), as well 
as coverage of the segments of M31 that fell into the 
ACS/WFC or WFC3/UVIS chip gaps in any single ex- 
posure. Smaller sub-pixel dithering enables well-sampled 
(Nyquist) images to be generated from a set of under- 
sampled single exposures. Unfortunately, all of these 
tasks were difficult to satisfy simultaneously, given the 
limited number of exposures possible, the need to operate 
WFC3 and ACS in parallel, and the geometric distortion 
present in all cameras. The distortion causes the angu- 
lar pixel-scale to vary with field position in the cameras, 
such that larger dither steps would produce a highly vari- 
able and generally non-optimal pattern of sub-pixel steps 
(once the integral portion of the step is subtracted) as a 
function of field location. Constructing Nyquist-sampled 
images thus requires keeping the total size of the dithers 
small, which conflicts with the large dithers needed to 
bridge the CCD-gaps. 

With these qualifications, our primary goal was to use 
dithers to produce Nyquist-sampled images in as many 
filters as possible. We were able to do this in F475W, 
F814W, and F160W. The F475W and F160W filters each 
require four dithers such that the fractional portions of 
the shifts map out a 2 x 2 pattern of 0.5 pixel steps. In 
contrast, only two dither positions are required in F814W 
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to achieve Nyquist-sampling, given its larger PSF. Al- 
though the use of a single diagonal dither to produce 
Nyquist sampling is less intuitive than traditional 2x2 
dither pattern, a single diagonal shift can produce frac- 
tional offsets of 0.5 pixels in both X and Y at the same 
time, allowing the two images to be interlaced to produce 
a new image with a rectilinear sampling-grid that is a 
factor y/2 finer than the native ACS/WFC sampling^] 
This approach is sufficient to ensure adequate sampling 
for F814W, but not F475W, which requires 2x finer sub- 
sampling. 

The dithers used to produce Nyquist-sampling in the 
three aforementioned filters are summarized in Table [2j 
The specific dithers were designed so that where possi- 
ble the correct fractional sub-sampling could be achieved 
in both ACS and WFC3 cameras simultaneously. The 
magnitude of the dithers is actually large enough so that 
small defects in the detectors can be stepped over, satis- 
fying a second goal of the dithering, but not so large that 
geometric distortion ruins the accuracy of the shifts. 

The design of the dither sequence was helped by the 
artful interplay of the location of the guard exposures 
and filter changes in the primary and parallel cameras, 
so that simultaneous sub-sampling dithers were not al- 
ways required in both instruments at the same time. 
Of the four sub-sampling dithers needed for F475W 
and F160W, one dither is "exact" for each filter, while 
three are done simultaneously in the two filters and pro- 
vide slightly non-optimal sub-pixel sampling. The sin- 
gle F814W sub-sampling dither occurs between a fil- 
ter change in WFC3/UVIS and is optimal for that fil- 
ter. In general, Nyquist-images can still be readily con- 
structed from dithers t hat stray sig nificantly from the op- 
timal sub-pixel steps ( |Lauer| [l999), a consideration that 
must be allowed for, even with an optimally-designed se- 
quence, since the HST Fine Guidance Sensor (FGS) has 
limited accuracy in executing the dithers. 

The dithers in the ACS F475W and F814W filters are 
not large enough to fill in the gap between the two CCDs 
in that camera. However, since every region that falls in 
the chip gap in one brick will be imaged again in another 
brick, completely filled coverage with the ACS can still 
be achieved. We thus have chosen to emphasize Nyquist- 
sampling in all ACS filters at the expense of depth in the 
small portion of the survey that falls into the chip gaps. 

Likewise, while the dithers in F160W are large enough 
to counter defects of one or two pixels in extent, the shifts 
are unfortunately not large enough to completely sam- 
ple over the "IR blobs" or the "death star" , which have 
characteristic sizes of 10-15 pixels and 51 pixels, cover- 
ing a total of ~700 pixels and ^2000 pixels, respectively 
(ISR WFC3 2010-06, Pirzkal 2010). We have found that 
the regions of the IR blobs can still produce adequate 
photometry, although with slightly higher uncertainties 
due to the 5-30% reduction in sensitivity in the area of 
the blobs. The only area then that completely lacks 
WFC3/IR data is the 0.2% of the field covered by the 
WFC3/IR "death star". Again, we consider the benefits 
of obtaining Nyquist-sampling over most of the survey 

23 A perfect analogue is the interlacing of black and white 
squares on a chess board. The black and white squares each form a 
regular grid with the sampling interval \/2x larger than the spacing 
of the interlaced squares of the chess board. 



in F160W to outweigh the sacrifice of a small portion of 
the area lost to unrepaired defects. 

In contrast, we could not obtain the four-point dither 
pattern required to construct Nyquist-sampled images 
in either WFC3/UVIS filter. However, we do not ex- 
pect this to be a significant scientific limitation, given 
that the WFC3/UVIS data are not crowding limited, and 
that any star detectable in WFC3/UVIS will also be de- 
tected in the Nyquist-sampled F475W images. Again, 
in contrast to the ACS/WFC, we did include a single 
large dither of ~ 37 pixels to cover the chip gap in both 
WFC3/UVIS filters, since we have essentially no later 
duplicate coverage of the fields with this camera. 

We also could not obtain the set of images needed to 
achieve Nyquist-sampling in the WFC3/IR F110W filter, 
due to the buffer-dump limits. However, we expect that 
every star that is detectable in the F110W image will also 
be detectable in either of the Nyquist-sampled F814W or 
F160W images. 

All dithers are tabulated in Table |2j A concern at the 
time of the Phase-II preparation of the program was that 
the accuracy of the dithers needed to achieve simultane- 
ous Nyquist-sampling in the F160W exposures in parallel 
with the F475W exposures would be difficult to realize 
in practice. In fact, due to an unfortunate y/2 error in 
the scale of the dithers, the first season of data (summer 
2010) was not properly Nyquist-sampled. This error was 
corrected in time for all subsequent observations, how- 
ever, and we find that the HST FGS is indeed able to 
perform the sub-pixel dithers to the level required. Fig- 
ure [4] shows the dithers realized in F475W during the 
Dec 2010 - Feb 2011 period of observation. This filter 
makes the strongest demands for dither accuracy of the 
three filters in which we will construct Nyquist-sampled 
images, but nearly all dithers fell within 0.15 pixel of the 
targeted offsets, which i s sufficient to achieve the desired 
degree of sub-sampling ( Lauer|[l999 ) . 

3.5. Tiling Strategy 

We use the basic exposure sequence above (two or- 
bits per pointing, one with WFC3/IR and one with 
WFC3/UVIS in primary, and ACS/WFC in parallel) as 
the foundation for a highly efficient tiling scheme. 

The primary tiling is based on WFC3/IR, which has 
the smallest FOV of all three imaging cameras. We ar- 
range the WFC3/IR pointings in a 3x6 grid, with <5" 
overlap among the WFC3/IR FOVs; due to distortion, 
the overlaps range from ~ 1.5" to ~ 5" for adjacent 
pointings. The WFC3/UVIS exposures use a small off- 
set of 1.655 pixels in X and 2.96 pixels in Y with respect 
to the standard UVIS-CENTER aperture, such that the 
same pointing center is maintained with respect to the 
IR-FIX aperture used to define the WFC3/IR pointings. 
The overlap between adjacent WFC3/UVIS exposures is 
~ 0.5 — 1', due to the camera's larger FOV. 

We cover the 3x6 "brick" with ACS by observing the 
two 3x3 halves of the brick with two orientations, 180° 
apart. Due to telescope constraints, observations in these 
two orientations are taken ~6 months apart, in two "sea- 
sons" . In the first observing season, the observations pro- 
duce a contiguous 3x3 half-brick of overlapping WFC3 
observations, and an adjacent 3x3 half-brick of overlap- 
ping ACS observations. Due to the large ACS FOV, 
most points within the survey region have ACS data 
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from two to four separate visits. In the second observing 
season, the telescope is rotated by 180°, and completes 
the 3x3 WFC3 pointings on the half-brick covered pre- 
viously with ACS. The switch in orientation leads the 
parallel ACS pointings to simultaneously cover the 3x3 
WFC3 pointings taken during the previous orientation. 

The resulting 3x6 brick has complete areal coverage in 
all three cameras, producing a contiguous area of 12 / x6.5 / 
in 36 orbits; exposure maps of the bricks are shown in 
Figure |5j Each pointing in the brick follows a consistent 
naming scheme, such that Field 1 is the pointing in the 
northeast corner of the brick, and Field 18 is the pointing 
in the southwest. The full naming scheme for each target 
position is of the form "M31-B##-F##-XXX", where "##" 
represents a two-digit number; the number after B indi- 
cates the brick number, the number after F indicates the 
field number within the brick, and the XXX indicates the 
camera on the WFC3 or ACS instruments (UVIS, WFC, 
IR). Note that parallels are named according to the area 
they overlap, rather than for the location of the primary; 
thus, ACS images labeled as "B01-F04" overlap WFC3 
images of the same name. 

The default orientation of the brick was set to allow 
the observations to be maximally schedulable in both 
180° orientations. The optimal orientations had large 
numbers of schedulable days in the summer (peaking in 
July) and the winter (peaking in January). While a sin- 
gle orientation for the whole survey would be preferable 
for producing minimal overlaps between adjacent bricks, 
upon consultation with STScI it was decided to adopt 
two default orientations for the survey. This choice al- 
lowed slightly longer scheduling windows for the obser- 
vations, reducing their impact on the scheduling of other 
programs. The 11 bricks closest to the center of the 
galaxy (Bricks 1-10, and Brick 12) are observed with 
ORIENT set to 69 (winter) or 249 (summer). The remain- 
der of the bricks in the outer galaxy are observed with 
ORIENT set to 54 (winter) or 234 (summer). Each brick 
is assigned to a unique proposal ID number (PID). The 
approximate corners of the NIR footprint of the bricks 
and their PIDs are given in Table [3] 

3.6. Spectroscopy 

The PHAT imaging described here is being augmented 
with extensive spectroscopy. Individual stars in M31 
are bright enough be observed spectroscopically, allow- 
ing us to measure the kinematics, metallicities, spectral 
types, and physical parameters of star clusters and mas- 
sive main sequence, asymptotic giant branch (AGB), and 
red giant branch (RGB) stars. 

The majority of PHAT spectroscopy to date has been 
carried out with DEIMOS on Keck (R = 6000, -6400- 



9100A), using a similar observ ing set-up as|Gilbert et al 



2006 



, |Kalirai et al. (2006a), and Guhathakurta et al. 
2006 ) . A total of 21 slitmasks covering much ol the HS'l 
footprint have already been observed, and have produced 
—5000 radial velocities accurate to a v = 5kms _1 down 
to I = 22 (~2 m fainter t han the TRGB ) ; prel iminary 
reductions are described in Dorman et al.| (2012). These 
spectra are being used primarily lor kinematic decompo- 
sition of different M31 sub-populations. Future analyses 
will constrain stellar metallicities and spectral types from 
the Ca triplet and other absorption features. These data 



also can be used to estimate ionized gas metallicities from 
Ha/[NII]/[SII] emission lines. 

Additional spectroscopic programs are underway, with 
the goals of: (1) measuring accurate HII region abun- 
dances using weak line methods; (2) measuring spectral 
types of hot stars using medium resolution spectroscopy 
at < 6000A; (3) measuring the masses and metallicities 
of stellar clusters. When possible, these various spectro- 
scopic programs also target serendipitous "interesting" 
sources from the PHAT survey, including PNe, AGN can- 
didates, and X-ray counterparts. We are also extending 
the star cluster an d planetary nebula survey of Caldwell 
et al.|fl2009[|201lD to PHAT targets, with more than 200 
new spectra obtained thus far. These various spectro- 
scopic programs will be described in future papers. 

4. DATA 

We now describe the present state of the PHAT ob- 
servations and pipeline for image processing, stellar pho- 
tometry, and astrometry. We also include assessments 
of the current data quality. Because this project is still 
actively acquiring data, we expect future data releases 
to incorporate on-going improvements in data process- 
ing and in the resulting photometric catalogs. We note 
planned upgrades and processing revisions throughout 
this section, for cases where substantial development 
work is already underway. 

4.1. Observations 

Figure [l]and Table [1] summarize the observations dis- 
cussed in this paper. We list bricks, their PIDs, and 
the range of dates over which data were taken for each 
orientation (i.e., the winter observations which populate 
WFC3 observations of the eastern halves of the bricks, 
and the summer observations which populate WFC3 ob- 
servations of the western halves). Observations for the 
program began July 12, 2010, and complete 6-filter cov- 
erage is available for Bricks 1, 9, 15, 17, and 21 as of 
Fall 2011. Brick 23 is nearly complete, but two fie lds ar e 
being re-observed due to guiding failures (Section 4.1.1). 
Our primary scheduling priority is to finish bricks, and 
we therefore anticipate that observations of the current 
"half-bricks" (2, 5, 8, 12, 14, 16, 18, and 22) will be com- 
pleted in the winter 2011/2012 observing season. Ob- 
servations are currently scheduled for the first halves of 
Bricks 4 and 6, as well. 

Figure [Hlshows false-color images (generated by Zolt 
Levay of !3TScI) of Bricks 1, 9, 15, and 21, which were 
the first four bricks completed in the program. Figure [7] 
shows a small portion of a false-color image for Brick 9, 
and reveals the rich information that can be seen at full 
HST resolution. 

4.1.1. Known Problems 

In a data set this large, there are bound to be glitches 
in the observations. Here we briefly point out known 
issues with some of the observations 

There were several problems affecting observations of 
Brick 23 (PID GO-12070). The Fine Guidance Sen- 
sors (FGS) failed to acquire guide stars during Visits 
03 and 13, making data from these observations unus- 
able. As a result, Brick 23 is currently missing high 
quality WFC3 observations for Fields 03 and 13, and 
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ACS observations for Field 16; the ACS parallels in Field 
06 appear to be usable. Data will be taken for these 
fields in the winter 2011/2012 observing season. In ad- 
dition, one WFC3/UVIS exposure for target M31-B23- 
F12-UVIS was not completely read off the HST recorder 
before being overwritten. As a result, half of one line is 
missing from the image. 

There are also currently missing observations in Brick 
16 (PID 12106). In one of our requested orientations, 
there is no guide star available for the pointing that cov- 
ers Field 17 with WFC3 and Field 14 with ACS. We 
therefore have to observe these two regions with the 
other orientation, and will have to point at each inde- 
pendently. The additional orbits needed for these obser- 
vations have been "borrowed" from the interface between 
the two brick orientations, where fields in Brick 13 have a 
large degree of overlap with Brick 15. We anticipate that 
data will be taken for these fields in the winter 2011/2012 
observing season. 

An FGS guide star lock failure also affected one ob- 
servation in Brick 1 (Visit 3, Orbit 1, of GO-12058 on 
2010-12-17). This target was subsequently re-observed 
on 2011-01-03. 

In addition to the above observational difficulties 
(which should be remedied in upcoming observations), 
visual inspection of the images shows that a small por- 
tion of them are slightly affected by a number of dif- 
ferent artifacts, most of which are well-known features 
of previous HST imaging observations (see HLA ISR 
2008-01 by M. Stankeiwics, S. Gonzaga, and B. Whit- 
more for ACS, and ISR WFC3 2011-16 by P. McCul- 
lough for WFC3/UVIS). Strongly over-exposed stars 
in ACS/WFC and WFC3/UVIS, for example, exhibit 
"bleeding-charge" tails, which are common to all im- 
ages obtained with CCDs. At the same time, bright 
stars also generate a number of scattered-light artifacts 
in ACS/WFC and WFC3/UVIS, such as weak out-of- 
focus pupil reflections or "ghosts" , which often manifest 
as pairs of highly-elliptical rings of light, somewhat re- 
sembling a pair of reading spectacles. Bright stars also 
generate more diffuse or complex "flares" of light, known 
as "dragon's breath" to users of these instruments. Both 
the dragon's breath and the "spectacles" can occur at 
angles of a few arcminutes away from the bright stars, 
which may not be even present in the image. These ar- 
tifacts should have minimal impact on our photometry, 
given that they affect only a trivial fraction of the survey 
area. Moreover, because our photometry uses local sky 
measurements, these additive effects usually have little 
impact on our photometry, beyond increased noise from 
the elevated background. 

While scattered light artifacts will be repeated in all 
images within a given exposure sequence, there are also 
transient events that typically affect only one image in 
the sequence. Surprisingly common events are long trails 
due to "space debris" passing through the field during an 
exposure. In some cases the trail may affect both CCDs 
in ACS/WFC or WFC3/UVIS. These are most likely 
small particles in orbit around the Earth — in many 
cases the object generating the trail is clearly out of fo- 
cus, but a point source will come into focus only when 
it is > 2 x 10 4 km away from the telescope. Again we 
have done no special processing to remove these trails or 
the spurious sources that they may produce; such pro- 



cessing is unnecessary for ACS, where all but the faintest 
trails will be flagged as cosmic rays, but additional image 
flagging could potentially be beneficial for WFC3/UVIS, 
where there are only 2 exposures per filter. The trails 
are typically extremely straight and uniform, and thus 
catalog-level contamination might be described with a 
fairly simple model. At the other end of the scale, we 
have also detected a few asteroids, which may be seen at 
different locations in several exposures within a sequence. 

The other transient scattered light artifact seen in our 
data is the illumination of ~l/3 of the WFC3/IR chan- 
nel when pointed near a bright Earth l imb. This effect 
has be en described and characterized in lDalcanton et a"D 
(2012), and is likewise thought to have very little impact 
on the photometry. 

Several exposures exhibit strong cosmic-ray "clusters," 
in which several dozen cosmic ray events form a tight 
cluster, superficially resembling a small star cluster. 
These events can be corrected with standard cosmic-ray 
repair algorithms; however, it may be worth noting that 
most of the data of the image affected may be lost over 
the extent of the cluster. 

Lastly, we note that errors in the transmission of the 
images to the ground in rare cases has resulted in the loss 
of a small number of pixels in a few images. The lost pix- 
els typically occur in a small contiguous segment of a few 
hundred pixels within a single row on the CCD imagers 
(we have not experienced any lost data with WFC3/IR). 
In all cases the lost data has been flagged in automatic 
reduction of the images, both in the data-quality image, 
and as anomalous values in the flattened image itself. 

4.2. Image Processing 

The first goal of our data processing is the production 
of a homogeneous catalog of 6-band photometry for all of 
the stars detected in our survey area. Obtaining such a 
catalog requires well-calibrated images that are cleaned 
of cosmic rays (CRs). 

The ACS camera has been well-studied and is well- 
calibrated in general. However, the camera has had is- 
sues with the bias level since the replacement of key 
electronic components during Servicing Mission 4. For 
WFC3, the calibration has been steadily evolving over 
the course of our survey. These issues have required 
us to go beyond the standard pipeline image process- 
ing when bias-correcting, flat- fielding, and flagging CRs 
in our imaging data. We now describe the image process- 
ing and CR rejection for the ACS and WFC3 cameras. 

4.2.1. ACS 

The ACS data were processed starting with the * . f It 
images from the HST archive. These images were flat- 
fielded by the OPUS pipeline. As our data come in ev- 
ery 6 months, multiple versions of the OPUS pipeline 
were applied to different portions of our data. Depend- 
ing on the observation date, versions 2010_2, 2010_3, and 
2011_le were used to generate the *.flt images. Since 
the ACS calibrations have been stable for several years, 
any changes are likely to be slight and are unlikely to 
affect the photometry. 

Due to issues with the ACS bias level (ISR ACS 2011- 
05, Grogin et al. 2011), the pipeline-processed images 
suffer from several problems. The first major problem is 
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that the ACS bias level shows a clear striping pattern, 
whose row-to-row fluctuations vary from image to image. 
These bias level variations are only apparent in images 
where the sky level was low; in such cases the row-to-row 
changes in the bias level have a larger amplitude than the 
typical sky noise. The second notable problem is that, 
in rows where the default pipeline bias subtraction was 
poor, the initial data quality extensions from STScI flag a 
high percentage of pixels as being affected by cosmic rays. 
We reset these CR-flagged pixels to allow us to determine 
the CR-affected pixels independently, after correcting for 
the bias problems. 

To correct for the bias striping, we used a destriping 
algorithm developed by Norman Grogin (csc2_kill .py, 
ISR ACS 2011-05). This algorithm attempts to correct 
images for striping, but is only helpful in cases where 
the bias striping can be accurately measured above the 
sky noise and where the bias striping corrections are not 
inadvertently correcting for real structure in the back- 
ground sky (say, in the bulge, where there are large gra- 
dients in the unresolved light). Therefore, after initially 
applying the de-striping algorithm, we then test the re- 
sults to be sure that the changes to the image are at the 
level appropriate for bias-stripe correction. To do so, we 
plot histograms of the difference between the final and 
corrected pixel values in one CCD column. If this distri- 
bution has an amplitude and width characteristic of the 
bias striping as described in the ISR, we replace the orig- 
inal *.flt image with the de-striped version. Otherwise 
the exposure is deemed to be unaffected by the striping, 
and the original exposure is unchanged. Roughly 50% of 
the exposures were kept as-is. 

After bias striping is corrected (if needed), we flag 
CRs as follows. For our longer exposures, the exposure 
times are similar enough to allow for reliable rejections 
based on high outliers from the median of the image 
stack. Therefore, all exposures longer than 50 seconds 
were put through the PyRAF routines tweakshif ts and 
multidrizzle, which flags CR-affected pixels using this 
median-image technique. We found that CR-rejection 
was far too aggressive when the short "guard" exposures 
were included in the image stack. We therefore handle 
CR rejection in the short exposures independently. All 
ACS exposures shorter than 50 sec onds are instead pu t 



through the IDL routine lacosmic (van Dokkum |2001| , 
which uses Laplacian edge-detection to flag pixels that 
show sharp edges associated with CRs, and which proved 
to be effective for removing obvious CRs from short sin- 
gle exposures. Once these steps are completed, the ACS 
data are ready to enter our photometry pipeline. 

4.2.2. WFC3 

While our ACS imaging required only a few minor 
changes to the pipeline-processed images, WFC3 is a 
sufficiently new instrument that its entire calibration has 
been in a state of flux during our survey, making homoge- 
neous calibration more of a challenge. Furthermore, our 
WFC3/UVIS data were obtained with just two exposures 
per filter, making CR- masking particularly difficult. 

The WFC3 data were calibrated starting with 
the raw images, because the UVIS and IR flats 
and distortion calibrations are continuing to evolve. 
All of our data were processed in PyRAF with 
calwf3 in STSDAS version 3.12. We used a sin- 



gle set of flat fields for all the images used is 
this release (f 110w_lpf It .f its, f 160w_lpf It . f its, 
f 275w_lpf It .f its, and f 336w_lpf It . f its, with header 
dates of 2008-12-09, 2008-12-09, 2009-04-24, and 
2009-03-31, respectively). For consistency, we also 
used a single set of distortion correction tables 
(210710_uvis_idc.f its, t20100519_ir_idc . f its, and 
u7nl8501j_idc.fits for WFC3/UVIS, WFC3/IR, and 
ACS, respectively). It is known, however, that the dis- 
tortion has changed as a function of time, for ACS at 
least (ISR ACS 2007-08); in future releases, we will be 
solving for time-dependent distortion corrections inte rnal 



to the PHAT data set, as described below in Section 4.6 
We anticipate updating all calibration images before each 
major rerun of the full data set. 

After flat-fielding, we flagged CRs in the calibrated 
WFC3 images. The calibrated WFC3/IR images are es- 
sentially free of CRs, as expected due to the many non- 
destructive reads taken during data collection. However, 
the WFC3/UVIS images, which contain only two expo- 
sures in each filter, were plagued by CRs. We attempt to 
mitigate the CR effects by running all WF C3/UVIS ex- 
posur es though the IDL routine lacosmic ( |van Dokkum 
2001), as was done for the short ACS guard exposures. 
We also process the images through the PyRAF routines 
tweakshif ts and multidrizzle using the minmed algo- 
rithm to flag CR-affected pixels. 

Unfortunately, even after these techniques are applied, 
the WFC3/UVIS data still contain some obvious CRs. 
More aggressive CR rejection was found to eliminate cen- 
tral pixels of stars, and thus we cannot pursue more 
aggressive CR rejection at the image-processing level. 
Instead of risking degradation of photometry for real 
stars, we cull CR artifacts from our photometric cata- 
logs, since they have poor fits to the PSF model and are 
anomalously "sharp". We are therefore able to cleanly 
remove CR-affected photometric measurements in our 
post-processing, as will be detailed in our description of 
photometry. Thus, while the residual CRs result in less 
attractive looking images, they do not affect our final 
photometry catalogs significantly. We plan to continue to 
experiment with ways to generate cleaner WFC3/UVIS 
images as the survey continues. 

After calibration, all WFC3 and ACS images are 
masked using the data quality (DQ) extensions of the 
*.flt images, which are created during image process- 
ing. The only modification was to accept data in the 
small I R bl obs (data quality flag 512), as described in 
Section [374} All images are then multiplied by the appro- 
priate pixel area map to allow for accurate photometry 
to be performed on the original calibrated and undrizzled 
images. 

4.3. Measuring Resolved Stellar Photometry 

We measure stellar photometry on stacks of the final 
calibrated and CR flagged images. The PHAT images are 
extremely crowded, making point spread function (PSF) 
fitting the only viable technique for producing accurate 
photometry. This technique requires a well-measured 
PSF model, calibration of that model against aperture 
photometry, and an algorithm for fitting the PSF and 
local background to all sources. 

All of our photometry was produced by t he software 
package DOLPHOT 1.2 QDorphin et al.|2002[ and many 
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unpublished updates). This software iteratively i denti- 
fies peaks and uses the PSF models from TinyTim (Krist 
1995; Hook et aT| 2QQ8) to simultaneously fit the sky and 
PSF to every peak within a stack of images, for mul- 
tiple filters taken with a single camera. Minor correc- 
tions for differences between the PSF model and the true 
PSF in each exposure are calculated by DOLPHOT using 
neighbor-subtracted bright stars in the field, primarily to 
account for changes in the telescope focus. Since each ex- 
posure is fit simultaneously, the original HST PSF is left 
intact to provide the highest possible quality PSF fits, 
making these corrections as sma ll as possible, typically 
of order 1-5% (see Section 4.3.1). 

The exact processing steps employed by DOLPHOT 
are described in detail in Dolphin (20001). To summarize, 



DOLPHOT uses the following steps to determine the star 
list and stellar photometry from a stack of aligned images 
for a single camera: 

• An iterative search for stars in the image stack. 

• Iterative PSF-fitting photometry of individual 
stars, using images with stellar neighbors sub- 
tracted. 

• PSF adjustment using bright, relatively isolated 
stars on neighbor-subtracted images. 

• A second round of iterative PSF-fitting photometry 
on neighbor-subtracted images, using the refined 
PSF. 

• Calculation and application of aperture corrections 
derived from bright, relatively isolated stars on 
neighbor-subtracted images. 

• Computation of individual VEGAMAG instrumen- 
tal magnitudes and uncertainties. 

• Application of CTE corrections to the ACS mag- 
nitudes; WFC3/UVIS CTE corrections will be in- 
cluded when they become available from STScI. 

• Computation of combined per-filter VEGAMAG 
magnitudes from all measurements that are unsat- 
urated and contain enough valid pixels near the 
core (i.e., do not fall on a cosmic ray or bad col- 
umn). 

Because DOLPHOT works on the uncombined image 
stack, precise alignment of the images is crucial for pro- 
ducing reliable PSF photometry. Any errors in image 
alignment would lead directly to systematic photometric 
errors, due to the PSF peak being misaligned with those 
of the stars. DOLPHOT solves for the relative align- 
ment of images by matching bright stars, starting from 
an initial estimate of the shifts provided by the user. The 
distribution of alignment values taken from the matched 
stars typically shows a clear central peak, with a width of 
order ~0.1 pixels for the star-to-star alignments within 
a single image. The resulting image-to-image alignment 
is therefore good to <C0.1 pixels. The PHAT processing 
pipeline flags cases where the distribution of alignment 
values lacks a central peak or has a width of a pixel or 
more; these rare cases are then subject to further by- 
hand intervention. 



We adopted DOLPHOT parameters (Table [§ that 
would give high completeness in highly-crowded regions 
and minimize systematics in each camera, and that also 
could be used for fields with a wide range of stellar den- 
sity (thus yielding homogeneous photometry from the 
bulge to the outer disk). We selected these initial param- 
eters ba sed on past experience with stellar population 
surveys (Dalcanton et al. 2009, 2012), and on a coarse 
sampling ol parameter space. The resulting parameters 
appear to have accomplished our goals, as they yield re- 
liable results throughout t he su rvey area according to 
our fake star tests (Section 4.7 below). However, given 



that we have identified a few issues affecting photometry 
at the few percent level (described below), we expect to 
optimize the processing parameters further as the sur- 
vey progresses. As such, we are currently experiment- 
ing with a larger grid of potential DOLPHOT parameter 
sets to quantitatively determine the optimal set for pre- 
cision and completeness. These parameters will likely 
be updated in our next data release, with details of any 
changes being described in subsequent papers and in sup- 
porting materials provided with the data distribution via 
the Multi-Mission Archive at Space Telescope (MAST). 

The key parameters for producing high quality pho- 
tometry have proved to be Rchi and Fit sky. Rchi is 
the radius inside of which the centroid of the PSF is fit. 
This value is set to be quite small (about one full-width 
at half maximum of the PSF) to avoid it being affected by 
poorly-subtracted neighbors. Subsequent star-matching 
between cameras and astrometric alignment also uses 
this same value of Rchi. Fit sky sets how the local back- 
ground level is measured; we have adopted Fitsky= 3, 
which forces simultaneous fits to the PSF and the back- 
ground level inside of the fitting radius. This strategy 
is optimal for very crowded fields, but is also computa- 
tionally expensive. For this first pass through our data 
we conservatively applied this algorithm for sky measure- 
ment. 

Once the PSF- fitted magnitudes of all stars have been 
measured, the values are corrected to a 0.5" aperture, 
equivalent to using an aperture correction measured us- 
ing nei ghbor- subtracted bright stars on the image (see 
Section [4.3.2] ). 

The AITS7WFC and WFC3/UVIS images both suf- 
fer from charge transfer efficiency (CTE) problems, such 
that not all charge is read out of a pixel, leaving resid- 
ual charge that is instead read out in subsequent rows. 
This effect leads to a small bleed trail of charge that 
is most severe for low electron levels and for rows near 
the beginning of the readout (i.e., near the chip gap for 
the ACS/WFC and WFC3/UVIS packages). The trail of 
charge is made up of flux that should have been included 
in the main body of the star, and thus photometry of 
stars affected by CTE will be systematically faint, and 
positions of the affected stars will be slightly shifted in 
the Indirection. 

Luckily, CTE is sufficiently well-behaved that it is pos- 
sible to derive reasonable corrections for its effects. The 
fluxes for ACS measurements are corrected for CTE us- 
ing the prescription of Chiaberge (ISR ACS 2009-01). 
Typical corrections range from to 0.2 magnitudes in 
Brick 01 where the background is high (hundreds of 
counts per pixel) to to 0.5 magnitudes in Brick 21 
where the background is low (tens of counts per pixel). 
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y-coordinates closer to the chip gap have larger correc- 
tions, as do fainter stars, which leads to the broad range 
in correction values. 

We have also examined the possibility of using the 
image-based CTE corrections of Anderson & Bedin (ISR 
ACS 2010-03), which correct images for CTE loss prior 
to running our pipeline. Unfortunately, the procedure 
in its current form has the characteristic that all peaks, 
including noise spikes in the background, are treated as 
stars. This has the result of magnifying background noise 
peaks and CRs above our detection threshold, creating 
enormous numbers of very faint false detections, most 
notably near the chip gap. Given this complication, we 
have chosen to retain the catalog-based CTE corrections 
for their superior noise properties at faint magnitudes. 

Currently, there are no CTE corrections available for 
WFC3/UVIS. Given that CTE is clearly evident in the 
WFC3/UVIS images (manifesting as obvious trails in Y 
for sources near the chip gap), our WFC3/UVIS photom- 
etry is undoubtedly biased to fainter fluxes for sources 
close to the chip gap. We unfortunately cannot con- 
strain the size of the effect directly with our data, since 
our tiling scheme is such that overlapping WFC3/UVIS 
observations leave stars at comparable Y chip positions 
(see Figure [5]). We will be incorporating WFC3/UVIS 
CTE corrections into our photometry as they become 
available. 

Finally, the fluxes are transformed to infinite aperture 
and converted to Vega magnitudes using the zeropoints 
dated 15-July-2008 in the ACS Users Handbook and ISRs 
WFC3 2009-30 and WFC3 2009-31 (J. Kalirai et al. 2009) 
for WFC3. These zeropoints will be updated should new 
values be published before our next run through the full 
data set. 

The photometry for this first pass through the PHAT 
survey data was performed on the data from each field 
and each camera separately. Therefore any depth that 
may be gained from overlapping data from neighbor- 
ing fields is not included in the current photometry. 
This added depth will be of only modest benefit in the 
crowding-limited portions of the survey, but will be sig- 
nificant in other portions. We anticipate that our second 
data release will contain simultaneous photometry on full 
image stacks, for all cameras. 

Before discussing t he pr oduction of the final photomet- 
ric catalogs (Section |4.4| ), we briefly discuss the ampli- 
tudes of the PSF and aperture corrections. 

4.3.1. Amplitude of PSF Corrections 

The TinyTim point-spread function models used for 
PSF photometry are unfortunately not perfect. For ACS, 
they are known to vary with time (see Figures 8 and 9 
in ISR ACS 2006-01 by J. Anderson a nd I. King) and 



ling) 

to have small errors with position (e.g., Jee et aL|2007 ) 
For WFC3, the errors are likely to be even larger, be 
cause no post-flight updates to the TinyTim models have 
been implemented at this time. 

Due to these temporal variations in the PSF, as well 
as possible systematic errors in the TinyTim PSFs them- 
selves, DOLPHOT computes a PSF adjustment image 
that is added to each PSF. This calculation is done mid- 
way through the iterative photometry process, once a 
reasonably final star list has been reached. The process 
involves analysis of residuals near bright, relatively iso- 



lated stars in images in which all detected stars have 
been subtracted. Although an entire difference image 
(actual minus TinyTim model) is calculated and applied 
to the photometry, only the value of this image in the 
central pixel is reported. A positive value indicates that 
the data are sharper than the model; a negative value 
indicates the reverse. 

Figure [9] shows the adopted PSF correction for the cen- 
tral pixel7~tor every image taken prior to Fall of 2011, with 
the exception of the short "guard" exposures in F475W 
and F814W. The amplitude of the typical PSF correc- 
tion is relatively small (<0.05 in the UV, and <0.01 in 
the optical and NIR), and has no strong dependence on 
the local stellar density (i.e., Brick number) or the time of 
observation (not shown). Table [5] lists the median PSF 
correction and the semi-interquartile width, as well as 
the number of stars used to calculate the PSF correction. 
The scatter in the measured corrections increases for the 
bluer cameras, most likely due to the smaller number of 
high signal-to- noise sources in the WFC3/UVIS images. 

The sign and amplitude of the central pixel PSF correc- 
tions are a first order indication of the errors in the PSF 
model. For the WFC3/UVIS data, the PSF corrections 
are substantial and are consistently negative, indicating 
the the current TinyTim model for the WFC3/UVIS PSF 
is slightly broader than the true data (i.e., the PSF cor- 
rection must put extra flux back into the central pixel). 
The PSF corrections for ACS and WFC3/IR are a fac- 
tor of three smaller, suggesting that the TinyTim models 
for these cameras are in better agreement with the data, 
at least within the very core of the PSF. The agreement 
outside the central pixel, however, is unconstrained by 
these initial tests. 

4.3.2. Amplitude of Aperture Corrections 

To calculate a star's magnitude, DOLPHOT adjusts 
the flux in the corrected TinyTim PSF to minimize resid- 
uals in the image stack. The model PSFs are scaled to 
unity within an aperture of 3" in radius, such that the 
internal magnitudes within DOLPHOT are calibrated 
as if the zero points were for 3" apertures. However, 
especially in crowded fields, we find it more practical 
to compute aperture corrections within a smaller aper- 
ture of 0.5" in radius, and then apply encircled en- 
ergy corrections to calibrate our photometry using the 
published zero points (which are for infinite aperture). 
DOLPHOT therefore corrects its internal PSF magni- 
tude to an aperture magnitude calculated within a 0.5" 
radius. DOLPHOT calculates these aperture corrections 
by identifying ~200 bright, reasonably isolated stars 
in each image, computing aperture photometry on the 
neighbor-subtracted residual image, and then measur- 
ing the differences between these aperture magnitudes 
and the stars' PSF magnitudes. These measurements 
are then combined in a weighted average with outlier 
rejection to derive a single aperture correction which is 
applied to the PSF magnitudes of all stars in the image, 
for a given filter. 

Figure [lO] shows the aperture corrections for ev- 
ery image taken prior to Fall of 2011, excluding the 
short "guard" exposures in F475W and F814W; Ta- 
ble [5] lists the median aperture correction and the semi- 
interquartile width for each filter, as well as the median 
number of stars used to calculate the aperture correc- 
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tion. The number of stars is typically high (~200) for the 
crowded optical and NIR data. The UV fields, however, 
frequently suffer from a paucity of bright sources, lead- 
ing to much smaller numbers of suitable stars and larger 
image-to-image scatter; this problem is most severe in 
Bricks 22 and 23 in the outer disk, and Brick 1 in the 
bulge. The Brick 1 bulge fields also suffer from reduced 
numbers of suitable stars in the optical and NIR, though 
not to the same degree as in the UV. The image-to-image 
scatter tends to be systematically larger for bluer filters, 
due to the lower number of high signal-to-noise stars. 

Table [5] shows that the mean aperture correction for a 
given filter is small (< 0.06 mag for all but F160W), and 
has a typical image-to-image scatter of less than 0.02 mag 
across the entire data set. However, inspection of Fig- 
ure [lO] indicates that the aperture corrections do have 
a modest systematic variation with brick number. We 
found no time-dependence in the amplitude of the aper- 
ture correction, suggesting that the observed variation 
with brick number is likely driven by the systematic de- 
crease in the stellar density and/or background level with 
increasing radius, since larger brick numbers and even 
brick numbers have larger projected radii. The effect is 
most noticeable in Brick 1, which has the largest range of 
stellar densities internal to the brick, due to the presence 
of the center of the bulge in Field 10. 



In contrast to what is seen in Figure 10 , the aperture 
correction should in principle be a constant that solely 
corrects for the difference in the encircled energy between 
0.5" and 3", assuming perfect PSF models and measure- 
ments. The fact that the aperture correction depends 
somewhat on the local level of crowding indicates that 
at least one of the key assumptions has been violated. In 
other words, the data in Figure [To| suggest that there is a 
slight bias in the measurement oFthe aperture correction 
calculated for bright stars, which may result from some 
combination of errors in the PSF model, background- 
dependent biases in the initial PSF magnitude, or sub- 
tle problems in constructing the empirical PSF used to 
measure the aperture correction. As such, it is not clear 
whether the observed crowding-dependence in the aper- 
ture correction is introducing systematic errors or cor- 
recting for them. 

At the present time, the ~0.02 mag scatter must be 
included in the error budget of the resulting photometry, 
with the understanding that there may be small system- 
atic offsets in photometry measured in different regions 
of M31, particularly in Brick 1 for the NIR. However, the 
systematic variations in the aperture correction are typ- 
ically smaller than empirical systematic errors derived 
from repeat measurements (shown below in Figure 18), 



suggesting the aperture corrections are not dominating 
the systematic error budget. They are also much smaller 
than the typical random uncertainties due to photon 
counting (Figure 13 , below). We are currently investi- 
gating the origin oTthis bias, and expect to have it and 
other small systematic uncertainties resolved before the 
second data release. 

4.4. Producing Photometric Catalogs 

The lowest level data products produced by 
DOLPHOT are "*.phot" files, which include the 
photometry for every detected source (both optimized 
for the entire image stack, as well as for each individual 



image), along with several metrics for the quality of 
each photometric measurement. These catalogs offer the 
greatest amount of information, but are frequently much 
larger than needed for most applications. We therefore 
generate additional smaller catalogs restricted only to 
stellar sources with good photometric quality. 

The metrics we have found most use- 
ful for generating high-completeness, 
low-contamination catalogs are measures of the 
"sharpness" and "crowding", which are quantified by 
DOLPHOT and returned as the parameters "sharp" 
and "crowd" for each star. Sharpness quantifies the 
degree to which an object is more o r less centrally 
peaked than the PSF (see eqn. 14 of |Dolphin| |2000| ) . 
Very centrally peaked sources (positive sharpness) are 
likely CRs, as they do not have the wings associated 
with true point sources. Very diffuse sources (negative 
sharpness) are likely background galaxies or unresolved 
groups of stars. Crowding quantifies the degree to 
which the photometry of the source may have been 
affected by neighboring sources. The value for crowd 
is determined by calculating the difference between the 
magnitude measured for the star before and after the 
subtraction of neighboring sources. Large numbers of 
bright neighboring sources therefore lead to a large value 
for the crowding parameter. 

Using these metrics, we cull lower-quality measure- 
ments from our catalogs to produce color-magnitude dia- 
grams showing well-defined features. The cuts were done 
at two levels for each camera, generating two catalogs 
with different restrictions on data quality. The first cuts 
the photometry at the most basic level, to ensure that 
all stars in the catalog have a high signal-to-noise de- 
tection in at least one filter. The resulting "st" cata- 
log contains only objects that either were detected with 
signal-to-noise S/N>4 in both filters or were detected in 
one filter at S/N>4 with sharp 2 <0.1; the additional 
sharpness cut is used for the one-filter detections to sup- 
press CRs. These cuts admit a very high percentage of 
all measurements but tend to remove any residual CRs 
from the catalogs. 

The second cuts are designed to produce sharp- 
featured CMDs with well- measured colors, by removing 
highly crowded stars, stars with poor PSF fits and stars 
with very low S/N in either filter. These cuts were made 
by keeping only measurements with S/N>4 in both fil- 
ters for a given camera, crowdi + crowd2 < 1.0 (0.48 
for WFC3/IR), and (sharpi + sharp 2 ) 2 < 0.075 (0.12 
for WFC3/IR). These cuts were determined by exper- 
imenting with a series of possible cuts, looking at the 
CMD of the rejected objects, and choosing the cuts that 
removed the largest number of stars without generating 
strong features in CMDs made from the rejected stars. 
These gst ("good star") catalogs contain only stars with 
the highest quality photometry, but suffer from more in- 
completeness in crowded regions, such as dense stellar 
clusters. They are also limited in depth by the require- 
ment that stars have good signal-to-noise in both filters 
for a given camera. In cases where the two filters have 
very different depths for typical stellar colors (such as for 
F275W and F336W), the gst catalogs will be missing 
many high quality detections in the deeper filter. Thus, 
when combining data from across cameras (say, when 
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generating a F336W - F475W CMD), it may be prefer- 
able to revert to the st catalogs to avoid this restriction, 
but then to apply joint crowding cuts on the detections 
from the different cameras. 

For PHAT data releases, the st and gst photometry 
catalogs are packaged into binary FITS tables, and dis- 
tributed via MAST. 

4.5. Artificial Star Tests 

We assess the completeness and uncertainty as a func- 
tion of color and magnitude of the st and gst catalogs 
using extensive artificial ("fake") star tests. We insert a 
series of 10 5 fake stars into each field with the appropri- 
ate PSF, for each camera independently. This strategy 
produces 5.4 million fake star tests per brick. Each fake 
star is added individually to the image stack (i.e., one 
artificial star at a time), and the photometry is rerun 
in a sub-region around the star to determine if the star 
was found, and if so, to determine what its measured 
magnitude was; this method is efficient, and does not 
run the risk of artificially changing the local crowding. 
The measured properties of the fake stars are then put 
through our quality cuts to determine which stars would 
be included in the final st and gst catalogs. 

The properties of the fake stars are chosen such that 
50% of the stars randomly sample the entire observed 
range of color and magnitude, and 50% sample the ac- 
tual color-magnitude distribution measured for the indi- 
vidual field. Spatially, the stars are distributed randomly 
throughout each field. While they therefore will broadly 
sample the photometric quality of the field, they will not 
necessarily be representative of the photometry in stellar 
clusters, which have higher than average local densities, 
but which cover only a tiny fraction of the area. 

Taking a coarse grid in color and magnitude, we calcu- 
late the completeness in each box of the grid and inter- 
polate to find the 50% completeness points. The limiting 
magnitude in the red filter is then the 50% completeness 
magnitude in the red band, measured at a color that is 
blue enough to ensure that non-detection in the red filter 
was the cause for low completeness (i.e., blueward of the 
corner at the faint end of the color-magnitude diagram) . 
The blue band limit is then the 50% completeness mag- 
nitude in the blue band, measured at a color that is red 
enough to ensure that non-detection in the blue filter was 
the cause for low completeness (redward of the corner at 
the faint end of the color-magnitude diagram). These 
50% completeness values are then tabulated as our mag- 
nitude limits in each field with the current processing 
version. 

4.5.1. Survey Depth 

Based on the 50% completeness limits, the faintest 
limiting magnitudes in the PHAT survey are found 
to be roughly 25.1, 24.8, 27.9, 27.1, 25.0, and 24.0, 
in the F275W, F336W, F475W, F814W, F110W, and 
F160W filters, respectively. For solar metallicity main se- 
quence stars, these depths correspond roughly to 2.6 M , 
2.8 M , 1.5 M , 1.5 M , 3.5 M , and 6M , respectively, 
for zero foreground extinction. However, these depths 
are a strong function of local stellar density for the opti- 
cal and NIR, as we now show. 

4.5.2. Crowding and Depth vs Radius 



For isolated stars, the depth of photometry is a con- 
stant set by the uncertainty in inferring a star's flux from 
a limited number of photons. At the surface density and 
distance of M31, however, one cannot assume that stars 
are isolated, or that their primary photometric uncer- 
tainties are due to photon counting. In M31's disk, the 
number of stars that fall above the limiting magnitude 
is sufficiently high that the PSFs of the individual stars 
will overlap. While crowded field photometry can com- 
pensate somewhat when detecting and measuring pho- 
tometry for partially overlapping stars, at some point 
the level of overlap is sufficiently large that faint stars 
are completely undetectable. The crowding also adds ad- 
ditional photometric uncertainties to the detected stars, 
since the reported photometry for a single star now de- 
pends on subtracting flux from neighboring stars as well. 
In these regimes, the limiting magnitude is no longer set 
by Poisson counting of electrons in the CCD, but in- 
stead is "crowding-limited". In crowding-limited data, 
the depth of the observations increases little with in- 
creased exposure time (although the photometry of the 
detected stars does improve), and only observations with 
higher angular resolution, or of lower surface density re- 
gions can significantly improve the limiting magnitudes. 
In M31, both the ACS and the WFC3/IR observations 
will be crowding limited over much of the disk, with 
the crowding limit being brighter in the higher surface 
brightness central regions. 

We can use the 50% completeness limits to determine 
the radial dependence of the survey depth of the survey. 
We assign a radius to each field by calculating the me- 
dian right ascension and declination of the detected stars, 
and then solving for the deprojected major axis length 
of an ellipse passing through that position (assuming an 
inclined disk with i = 70°, a position angle of 43°, and a 
center at a(2000) = 10.68473° and £(2000) = 41.26905°; 
we discuss the c hoice of inclination and position angle 
further in Section |5^2| below) . We also calculate the mean 
number surface density of stars with S/N>4 (i.e., stars 
per arcsecond 2 ) in the st and the culled, high-quality 
gst catalogs for each field. 

We plot the radial dependence of the 50% complete- 
ness magnitude limit in Figure [TT] As expected, we see 
a strong radial dependence in the completeness limit, in 
spite of all fields having nearly identical exposure times. 
In both the optical and the NIR, there is a brighter mag- 
nitude limit in regions with higher stellar densities. The 
effective magnitude limit of the ACS filters is more than 4 
magnitudes brighter in the bulge than in the diffuse outer 
disk. The variation in the WFC3/IR channel is even 
larger, being ~5.5 magnitudes deeper in the outer disk 
than in the bulge; given the larger pixels of the WFC3/IR 
camera, it is not surprising that the NIR observations are 
the most severely affected by crowding. 

In contrast, there is little radial variation in the 50% 
completeness limit in the WFC3/UVIS data, with only 
the slightest hint of an upturn in the central bulge field. 
The limiting magnitude of the WFC3/UVIS data are 
therefore set by photon-counting, not crowding. Longer 
WFC3/UVIS observations could therefore be expected 
to yield significantly deeper CMDs, unlike for the PHAT 
WFC3/IR and ACS/WFC CMDs, which are the deepest 
that are currently possible at these radii in M31. We note 
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that the observed radial variation a grees extremely we ll 
with predictions initially made using Olsen et al. (2003)'s 
crowding simulations (modified for the PHAT observing 
parameters) in advance of the start of survey operations. 

We can see the impact of stellar density on photomet- 
ric depth more directly in Figure [l2j where we plot the 
correlation between the 50% completeness limit and the 
stellar density of all stars detected with S/N>4 in a given 
filter (open circles, from the st catalogs), and for the sub- 
set of stars with reliable photometry and S/N>4 in both 
of the filters observed with a single camera (solid circles, 
from the gst catalogs). At high stellar densities, the pho- 
tometric quality cuts significantly reduce the numbers of 
stars in the gst catalog compared to the less stringent 
st catalog, by roughly a factor of three. These cuts have 
a particularly large impact on the F336W filter, because 
of the requirement that the gst catalog also have S/N>4 
for the shallower F275W observations. 

For photon-limited observations, we expect the limit- 
ing magnitude in Figure [12] to be constant for a fixed 
exposure time, and to have no dependence on the stel- 
lar surface density. We see this expected behavior in 
the WFC3/UVIS filters, where there is no correlation 
between the completeness limit and the stellar density. 

The relationship between stellar density and the mag- 
nitude limit is far more complicated for the ACS data, 
which is crowding-limited over most of the disk. At 
the lowest stellar densities shown, there is no obvi- 
ous relationship between the completeness limit and 
the stellar surface density, indicating that some of the 
ACS data has reached a low enough surface density 
(< 6 stars arcsecond -2 ) that the photometry becomes 
photon-limited. These limiting magnitudes are roughly 
27.9 in F475W and 27.1 in F814W for the ACS expo- 
sures taken at a single pointing^] At higher stellar sur- 
face densities, the limiting magnitude of the ACS data is 
systematically brighter, as the crowding begins to affect 
the ability to detect faint stars (the st points), or to pro- 
duce reliable photometry (the gst points). Instead, the 
completeness limit becomes pinned to whichever depth 
produces ~ 14 stars arcsecond -2 with reliable photome- 
try (i.e., in regions with intrinsically high stellar surface 
densities, the density of bright stars is sufficiently high 
that the maximum surface density can be reached with 
only the brightest stars in the luminosity function, lead- 
ing to a bright limiting magnitude). 

The NIR data in Figure [12] is similar to the opti- 
cal ACS data in many respects. However, it lacks the 
turnover at low stellar surface densities to a single 50% 
completeness magnitude, indicating that the WFC3/IR 
data is crowding-limited throughout the disk, as was ex- 
pected from our initial simulations. The surface den- 
sity of reliably measured stars is essentially constant at 
~ 2 stars arcsecond -2 throughout the galaxy. The small 
ripples in the stellar surface density at ~24.3mag in 
F110W and -23.2 mag in F160W are due to the in- 
creasing numbers of stars at the magnitude of the well- 
populated red clump, and the bright ripples at 19 and 
20th magnitude (in F110W and F160W, respectively) 
are due the increasing numbers of stars at the NIR tip of 
the red giant branch (TRGB). Hints of a bump caused 

24 The overlapping ACS observations from adjacent pointings 
will increase the depths in these photon-limited regions. 



by the red clump is also visible in the ACS data (at 
~24.8mag in F814W), but it is less pronounced, because 
the red clump is somewhat blurred by dust in the opti- 
cal, making it a less pronounced feature in the luminosity 
function. 



Figure 12 also shows the substantial penalty that re- 
sults from the larger WFC3/IR pixel scale. At the crowd- 
ing limit, the ACS camera produces reliable photome- 
try for ~7 times as many stars per arcsecond 2 than the 
WFC3/IR channel. Deeper NIR data in M31 will only 
be possible with a telescope+camera system with smaller 
pixels and higher angular resolution. 

4.6. Astrometry 

Producing accurate 6-filter photometry for crowded- 
field photometry requires exquisite sub-pixel alignment 
between overlapping images. We have made care- 
ful astrometric solutions to place our PHAT data 
on a universal reference frame with an absolute er- 
ror of 100 milli-arcseconds (dominated by errors in the 
ground-based reference catalog) and relative errors of 
10 milli-arcseconds between stars observed in the three 
cameras. The astrometric alignment is a multi-stage 
process in which individual exposures within a field are 
aligned to each other, overlapping regions within bricks 
are used to tie the fields together, and finally bricks are 
aligned to a ground-based absolute astrometric reference 
frame. We describe each of these steps below. 

4.6.1. Alignment of individual exposures taken with one 
camera during one visit: 

As discussed briefly above, the individual exposures in 
an image stack are aligned within the DOLPHOT pro- 
cessing module. DOLPHOT identifies bright stars within 
each exposure, and then uses the geometric distortion 
correction to calculate relative shifts to apply to each 
exposure, based upon initial guesses generated from the 
astrometry in the header. Typically hundreds to thou- 
sands of stars are used in the alignment, with ~300, 
525, 850, 800, 6000, 6250 average alignment stars for the 
F275W, F336W, F475W, F814W, F110W, and F160W 
filters, respectively; the number of alignment stars used 
for the NIR filters varies dramatically with radius, from 
~2000 in the outer disk to 20,000 in the inner bulge. 
The resulting corrections to the (relative) header-based 
astrometry are typically of the order of 1-2 pixels or less. 
The only exceptions found in the data so far are for the 
ACS images in Brick 21, Field 2, which have relative 
shifts that differ substantially from what would be in- 
ferred from the headers; data quality in this field is oth- 
erwise unexceptional. For WFC3/UVIS and ACS/WFC, 
each chip is aligned independently. The RMS of the 
shifts for pairs of alignment stars are typically ~0.2 pix- 
els for WFC3/UVIS and ACS/WFC, and 0.1 pixels for 
WFC3/IR. 

4.6.2. Aligning all images for a single camera: 

The photometric catalogs initially generated by 
DOLPHOT contain distortion-corrected astrometry, tied 
to the astrometry of a reference frame generated with 
MULTIDRIZZLE. These reference frames, in general, are 
translated and rotated with respect to a global reference 
frame (due to errors in guide star catalog positions) and 
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also have shear, scale, and distortion (due to residuals 
from the pipeline distortion corrections). We aim to re- 
move these errors by aligning the individual catalogs with 
each other and with an astrometric reference catalog tied 
to a global reference frame. For aligning the astrometry 
we use the gst catalogs for ACS/WFC and WFC3/IR, 
since they contain fewer spurious stars, and the st cata- 
logs for WFC3/UVIS, with cuts on the F336W detections 
only (S/N>4 and crowd <0.25) to avoid eliminating the 
large number of valid F336W-detected stars that were 
undetected in the shallow F275W imaging. We split the 
ACS/WFC catalogs along the CCD gap and handle the 
two CCDs separately. 

We begin by coarsely aligning the individual catalogs 
in celestial coordinates, allowing only translational cor- 
rections. For the large ACS/WFC catalogs, we cut to 
bright (F475W < 24) stars for speed. Given a set of spa- 
tially overlapping catalogs, we search each pair of cata- 
logs for pairs of stars (one from each catalog) within 2" 
of each other, with this limiting separation set by the 
expected astrometric uncertainty of the cameras relative 
to the guide stars. This large matching radius is nec- 
essary to handle the typical offsets we find in our data. 
If the affine (rotation, shear and scale) and distortion 
corrections are small, then the distribution of positional 
differences between the pairs of stars will include a peak 
(due to true matches) and a background (due to false 
matches). We also observe a ring of lower match density 
around the peak due to DOLPHOT's closeness-of-sources 
limit. The position of the peak in (ARA, ADec) space 
is the mean positional offset between the two catalogs, 
and the width of the peak is due to measurement errors 
plus affine and distortion errors. We fit for (ARA, ADec) 
using a Gaussian peak plus flat background mod el using 
the Expectation Maximization algorithm (EM; |Demp-| 
ster et al. 1977). Once we have measured the offsets 
between each pair of overlapping catalogs, we solve a 
least-squares equation to find a consistent set of transla- 
tional corrections that will simultaneously minimize the 
distance between matched pairs of stars for the entire set 
of catalogs. 

4.6.3. Tying photometric catalogs to a global astrometric 

frame: 

After applying the initial translational corrections de- 
scribed above, the contiguous catalogs for each camera 
are roughly aligned internally. Next, we align the cam- 
eras with each other and with a global reference frame. 
This step is required for eventual comparisons with other 
mult i- wavelength data sets. 

We have defined a dense network of global astrometric 
standards using archival i' data from MegaCam CFHT 
(courtesy of Jean-Charles Cuillandre). The CFHT pro- 
cessed images were subjected to iterative PSF photom- 
etry. The photometry in the resulting astrometric cat- 
alog extends down to %' = 21. 5 and contains 1.7 m illion 
sources. For comparison, the |Massey et al. ( |2QQ6| ) Lo- 
cal Group Survey (LGS) catalog contains one-seventh as 
many sources over the same region. The positions of the 
brighter, non-saturate d stars were matched to 2MASS 
20061)- 
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-which itself is tied to Tycho-2 
-to produce a well-constrained astro- 
metric solution tor the entire region. The dispersion for 
matches between CFHT sources with i' < 18 and 2MASS 



is about 0.2 arcseconds. 

We begin the global alignment by roughly aligning the 
ACS/WFC catalogs to the CFHT-based reference cat- 
alog. We first merge the ACS/WFC catalogs, remov- 
ing one star from each pair within 30 milli-arcseconds of 
each other, and cutting to F814W < 21.5 to match the 
depth of the reference catalog. Then, as before, we find 
stars in the ACS/WFC and CFHT catalogs within 1" of 
each other, fit a Gaussian peak plus background model, 
and use the peak weights to fit a least-squares affine cor- 
rection (translation, rotation, scale, and shear) to the 
ACS/WFC catalogs. Similarly, we roughly align the 
WFC3/UVIS and WFC3/IR catalogs to the ACS/WFC 
catalogs. 

Finally, we repeat the alignment process, this time 
allowing each individual catalog an affine correction to 
bring it into alignment with its peers and with the ref- 
erence catalog. We do this, as before, by finding nearby 
pairs (both between pairs of catalogs, and between the 
catalogs and the reference), fitting a Gaussian peak plus 
background model to properly weight the matches, and 
computing a least-squares set of affine corrections j^j 

We characterize the accuracy of the photometry as 
the width of the Gaussian peak fit to the distribution 
of all offsets between pairs of matched stars. The final 
alignment of the ACS/WFC catalog to the CFHT refer- 
ence frame has a 1-cr accuracy of ~ 45 milli-arcseconds 
per coordinate in the outer disk, increasing to ~ 
60 milli-arcseconds in the bulge. The internal alignments 
within the WFC3/UVIS, ACS/WFC, and WFC3/IR 
catalogs have 1-cr accuracies of ~ 8.5 milli-arcseconds, 
~ 6 milli-arcseconds, and ~ 6 milli-arcseconds respec- 
tively. The relative alignments of the WFC3/UVIS and 
WFC3/IR catalogs to ACS/WFC have 1-cr accuracies of 
~ 6 milli-arcseconds and ~ 3.7 milli-arcseconds respec- 
tively. 

Even with the inclusion of affine corrections, we see 
clear structured residuals of up to ~ 5 milli-arcseconds 
when matching different cameras to each other. These 
residuals have a consistent structure from field-to-field 
within a brick, and are almost certainly caused by errors 
in the geometric distortion correction. In future releases, 
we expect to refine the above procedure to include cor- 
rections to this geometric distortion. 

4.6.4. Merging photometric catalogs for all cameras: 

For some applications it is helpful to have catalogs of 
matched multi-camera (6-filter) photometry. This step 
involves merging measurements made with three dif- 
ferent cameras, with multiple overlapping pointings, in 
highly crowded fields. As a result, this step can only 
be done after all catalogs have been tied together as- 
trometrically to a high degree of accuracy. We gen- 
erate merged, six-filter catalogs as follows. We first 
remove duplicates from catalogs of overlapping fields, 
for each camera. First, we identify and group any 
ACS/WFC stars that appear within 50 milli-arcseconds 
of each other in multiple catalogs. We then replace 

25 We note that we have more than ample signal to perform 
these fits, in spite of the large number of parameters needed for 
the field-by-field affine transformations. The least-squares problem 
when aligning a set of 18 ACS/WFC catalogs to the CFHT-based 
reference typically has ~ 10 7 matches, but only 216 parameters (6 
per catalog). 
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them with a single ACS measurement assigned to the 
mean of the measured positions and magnitudes in the 
overlapping catalogs. We do the same for WFC3/UVIS 
and WFC3/IR, resulting in a single brick- wide cata- 
log for each camera, with duplicates removed. We 
then generate the multi-camera photometric catalog by 
matching the merged ACS/WFC stars with the merged 
WFC3/UVIS and WFC3/IR stars, using a matching ra- 
dius of 50 milli-arcseconds. Stars are not required to be 
matched in multiple cameras, such that a source with 
no counterparts within 50 milli-arcseconds will still be 
included in the merged catalog. 

We note that deblending sources near the detection 
limit in highly crowded fields does not always produce 
a unique result, leading to slight variations in the po- 
sitions and assigned magnitudes of deblended stars. As 
such, source matching becomes less reliable at faint mag- 
nitudes. In future releases, we will be doing source iden- 
tification and photometry simultaneously across all over- 
lapping images in all three cameras, which should elimi- 
nate this problem. 

4.6.5. Propagating astrometry back into the photometric 

catalogs: 

The result of our astrometric alignment process is an 
affine transformation for each camera in each field. We 
apply these affine transformations to the FITS World Co- 
ordinate System (WCS) headers of the input *.flt im- 
ages, along with the small relative shifts that DOLPHOT 
calculated during alignment of the image stacks. We 
modify the WCS CRVAL headers to apply the shifts, and 
the CD matrix elements to apply the remaining affine 
terms. We also propagate the revised astrometry back 
to the individual photometric catalogs. 

4.7. Data Quality 

We have taken a number of different steps to character- 
ize the current quality of the photometry. In the follow- 
ing sections we discuss the raw p hotom etric uncertainties 
reported by DOLPHOT (Section 4.7.1 ), the uncertainties 
and biases du e to c rowding, as assessed by artificial star 
tests (Section 4.7.2[ ) , and uncertainties due to systematic 
errors, as measured from repeatability tests in overlap- 
ping images (Section 4.7.3). All results in this section are 
for individual fields in the current data release, which do 
not take advantage of additional depth due to overlap- 
ping adjacent fields, or upcoming improvements in the 
flat fields, distortion corrections, or PSF models. We 
therefore expect that some of these uncertainties will be 
reduced in upcoming versions of the PHAT photometry. 

4.7.1. DOLPHOT -reported Uncertainties 

DOLPHOT reports photometric uncertainties for each 
star as part of its standard output. These uncertainties 
are based on the Poisson variation in the flux from the 
star, from the sky, and from neighboring stars whose flux 
has been subtracted. The value of the uncertainty (in 
magnitudes) can be found in the MAG1_ERR and MAG2_ERR 
entries of the gst and st photometric catalogs. 

In Figure [l3| we plot the median photometric uncer- 
tainty as a function of magnitude, for all 6 filters, in Field 
9 of Bricks 1, 9, 15, and 21 (upper left to lower right). 
The median uncertainty is given in magnitudes, and was 



calculated for stars with S/N>4 in the named filter, using 
data from the st catalogs. The median uncertainty and 
magnitude was calculated for groups of 25, 400, and 200 
stars, for the WFC3/UVIS, ACS/WFC, and WFC3/IR 
cameras, respectively. We have not included stars that 
only appear in a subset of the images (i.e., those that are 
near the edges of an image or near a chip gap); such stars 
have much higher reported photometric uncertainties, as 
expected, but are not representative of the bulk of the 
photometric catalog. We have also excluded stars that 
touch the edge of an image, or that have saturated pixels 
(i.e., FLAOO). 
The median photometric uncertainty in Figure 13 cor- 



relates strongly with magnitude, as expected. However, 
the slope of the correlation is different for each camera, 
due largely to the changing importance of the sky bright- 
ness. In the UV, the uncertainty is dominated by the 
Poisson variation in the flux from the stars themselves, 
since the PHAT WFC3/UVIS images are essentially un- 
crowded and the background of unresolved flux is quite 
dark. In the optical and NIR, the situation is quite dif- 
ferent, due to a high level of crowding and brighter back- 
grounds. The correlation between photometric uncer- 
tainty and magnitude therefore follows a different slope 
than in the UV. Moreover, the exact relation changes 
with radius as well, due to the darker background and 
lower stellar density at large radii. 

We note that the photometric uncertainties shown in 
Figure 13 are not always the dominant source of error. 
For faint stars, flux from unresolved stars below the de- 
tection limit can significantly bias the resulting magni- 
tudes. For bright stars, systematic uncertainties in the 
flat field, CTE corrections, or PSF model can dominate 
over the formall y sma ll photometric error. As we show 
below in Section [4. 7. 3[ these systematic errors have typ- 
ical amplitudes of 0.02-0.04 magnitudes for bright stars; 
we include a horizontal line in Figure 13 to indicate the 
approximate magnitude range where systematic uncer- 
tainties will become important. 

4.7.2. Crowding Uncertainties and Biases from Artificial 

Stars 

The accuracy of stellar photometry can be strongly af- 
fected by crowding in regions of high stellar density. In 
such cases stars that are too faint to be detected individu- 
ally can blend with other stars to rise above the detection 
limit. Thus, many of the faint stars in the photometric 
catalog have been biased to brighter magnitudes by their 
neighbors. High crowding also makes it difficult to ac- 
curately measure the local background, because all of 
the pixels outside a star's photometric aperture are filled 
with other stars. Finally, crowding tends to increase the 
photometric uncertainty, because the measured flux de- 
pends on the accuracy with which the flux from all sur- 
rounding stars can be subtracted from the image. 

The only way to accurately assess the effect of crowd- 
ing is through e xtens ive artificial stars tests (described 
above in Section 4.5). These tests allow us to measure 



the difference between the true and recovered magnitude, 
as a function of the star's color, magnitude, and position. 

In Figures [T4] - [T6l we show a series of three-panel plots 
containing the cumulative distributions of the magni- 
tude difference (left panels), the fractional flux difference 
scaled by the flux uncertainty (center panels), and the 
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color difference (right panels) as a function of magnitude 
(where redder colors and thicker lines indicate fainter 
stars), for the UV (Figure [l4|) , optical (Figure (15]), and 
NIR (Figure 16), using Fiela 9 of Bricks 1, 9, 15, and 
21 (top to bottom, respectively). We only include data 
for Brick 1 in the UV analysis; the WFC3/UVIS data is 
uncrowded at all radii outside the very central regions of 
the bulge, and thus the artificial star tests yield identical 
results at all radii. All quantities are given as the mea- 
sured value minus the true value for the inserted star. 
Photometry only includes stars from the gst catalogs. 

Figure [14] shows that the effect of crowding is mini- 
mal for the F275W and F336W filters. The differences 
between the input and output magnitudes and/or colors 
increase for fainter stars (left and right panels, respec- 
tively). However, if one scales the flux difference by the 
photometric uncertainty reported by DOLPHOT (center 
panel), it is clear that the distribution of flux differences 
is essentially Gaussian, with a width set by the photon- 
counting uncertainties alone. The only sign of crowding 
effects is an extremely slight tail towards brighter fluxes 
at faint magnitudes in the deeper F336W filter (thick red 
line in the central panel of the bottom row). 

The effects of crowding are more severe in the opti- 
cal, as can been seen in Figure 15 Unlike in the UV, the 
distributions of flux differences (center panels) do not fol- 
low the Gaussian distribution expected for uncertainties 
dominated by photon- counting (with the possible excep- 
tion of the brightest stars in the least crowded fields). 
The distributions are instead skewed, such that stars are 
preferentially detected as brighter than their true magni- 
tude. The skewing yields a significantly increased proba- 
bility that a star will be reported as being brighter than 
expected for photon-counting uncertainties (i.e., compare 
the fraction of stars that are more than 5<j brighter to 
either the Gaussian expectation, or to the fraction that 
are 5cr fainter). However, the median magnitude differ- 
ence is not highly biased outside of the bulge, and the 
median flux errors are typically biased by less than the 
expected photometric error. The effects of crowding are 
noticeably worse in the inner regions of the galaxy, where 
the surface density is highest. Crowding also has more 
of an effect in the F814W filter, which has a larger PSF, 
and which is deeper than F475W, and thus has a higher 
density of resolved stars. 



The crowding errors are largest in the NIR (Figure 16), 
which suffers from worse resolution than either ACSor 
WFC3/UVIS, due to bigger pixels and an intrinsically 
broader PSF at long wavelengths. As a result, the effects 
of crowding are severe at all radii, even the outer disk. 
These effects can be clearly seen in the distributions of 
magnitude errors, which are dramatically wider than ex- 
pected from photon counting along (middle panels), and 
which show clear skewing (left panels). 

Surprisingly, the NIR photometry is biased towards 
fainter magnitudes for the median star. The bias is slight 
in F110W, but noticeable in F160W, where the crowd- 
ing is most severe. This difference between the two fil- 
ters then produces a slight blueward shift in the median 
F110W-F160W color (<0.05 mag in all but the faintest 
bin). In all cases, the systematic bias in the median mag- 
nitude or color is significantly smaller than the width of 
the error distribution, such that the error properties will 
be dominated by random rather than systematic errors 



in most applications. The systematic bias is also of the 
same order as both the in trinsic errors in repeat measure- 
ments (see Section 4.7.3| below) and the likely uncertain- 
ties in the current zero points. As such, while the biases 
in Figure [16| are undesirable, they are unlikely to be the 
dominant source of systematic errors, except perhaps in 
the densest regions in the bulge. Unfortunately, we do 
not currently have a convincing explanation for this off- 
set at this point. We are continuing to investigate, and 
expect to have the issue resolved for the second data re- 
lease. 

4.7.3. Repeatability and Systematic Uncertainties 

There are a number of position-dependent uncertain- 
ties that are not reflected in the DOLPHOT-reported 
uncertainties or the artificial star tests, but that will con- 
tribute to magnitude differences when the same star is 
measured in different parts of a single chip. These include 
errors in the flat-field, the position-dependent TinyTim 
model for the PSF, the CTE correction, and the image- 
to-image DOLPHOT aperture corrections. These will 
manifest in the distribution of magnitude differences as 
a constant systematic error that adds linearly with the 
random photometric error. We can test for these effects 
via the repeatability of our photometry when comparing 
the magnitudes of individual stars in overlapping obser- 
vations. 

We use the PHAT astrometry to identify all pairs of 
matched stars from the gst catalogs whose positional off- 
sets give them a better than 75% chance of being a true 
match, based on the observed distribution of astrometric 
offsets. These matched pairs sample the entire area of 
the ACS/WFC, due to the large amount of overlap be- 
tween adjacent fields. In the WFC3/UVIS camera, the 
matched pairs cover more than 50% of the chips, but do 
not sample a square "hole" in the center of the field of 
view. In the WFC3/IR channel, the matched pairs only 
sample the very edges of the field of view (~ 3-5"), due 
to their limited overlap (see Figure |5b. Thus, one should 
be aware that only the results for ACS /WFC can be seen 
as truly representative of the entire field of view, whereas 
for the WFC3/IR channel, we are only probing the ex- 
treme edges of the chip, where many of the systematic 
effects are likely to be at their worst. 

In Figure [17] we plot the cumulative distribution of ob- 
served magnitude differences between repeated observa- 
tions of stars within limited magnitude ranges (i.e., bins 
of 0.5 magnitude, with bluer lines indicating brighter 
magnitude ranges). We also plot the la widths of the 
distributions as a function of the median magnitude in 
each bin (black line), derived from the observed inter- 
quartile range of the observed distribution, assuming a 
Gaussian distribution. These widths are compared with 
the expected width of the distribution, based solely on 
photon counting statistics (blue line). The quadrature 
difference between the observed and expected width (red 
line) is then a measure of the systematic error (although 
an imperfect one, since systematic errors are unlikely to 
be strictly Gaussian). These measurements are shown 
for Brick 15, but the results are representative of what 
we see in other disk fields with very different crowding 
levels. 

Figure [17] shows several clear trends. First, the dis- 
tribution of magnitude differences becomes systemati- 
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cally wider at faint magnitudes. However, comparing 
the observed to the expected width (black versus blue 
line) shows that the majority of this trend is driven by 
the larger variation in photon counting for small fluxes. 
Second, at bright magnitudes, the observed magnitude 
variations are much larger than expected from the for- 
mal photometric errors, by factors of 10 in some cases. 
This difference is a clear signal of systematic errors (red 
line). However, the systematic errors are small (< 0.05 
mag) at bright magnitudes, where it is easiest to assess 
their amplitudes. Finally, photometric errors begin to 
dominate at faint magnitudes in most cases. The magni- 
tude where photometric magnitudes begin to dominate 
depends on the camera, filter, and degree of crowding. 
For example, the WFC3/IR observations are crowding- 
limited, and thus can only detect stars much brighter 
than th e true photometric limit of the observations (see 
Section [4X2] ). Thus, the WFC3/IR detected stars tend 
to have lower average photometric uncertainties, making 
the systematic uncertainties proportionally more impor- 
tant. 

We have calculated the width of the observed distri- 
bution of color differences in a similar fashion. Because 
spatial-dependent biases in flat fields and PSF models 
are often correlated between filters, systematic biases in 
colors may in fact be smaller than those for magnitudes 
alone p*] In Figure 18 we compare the inferred level of 
systematic bias in color (black line) to those in the mag- 
nitudes of individual filters (blue and red lines, for the 
bluer and redder filters used for an individual camera). 
The systematic errors in colors are indeed smaller than 
those for individual magnitudes, by up to a factor of two 
in the NIR. They have typical amplitudes of only 0.02 
mag, rising to ~0.04 mag in the faintest bins. They also 
appear to have a weaker dependence on the brightness 
of a star, and are thus more constant with magnitude. 

We can use stars that are measured in multiple images 
to assess whether the systematic biases in Figure [l8| have 
spatial structure, as one would expect for flat-fielding er- 
rors or spatially-dependent errors in the PSF model. For 
ACS, we have the benefit of highly overlapping images, 
such that an individual star appears in multiple quad- 
rants of the ACS chips, and such that every position on 
each chip has overlapping data. We use these overlaps to 
assess whether there are locations in the ACS FOV where 
stars are consistently measured to be brighter or fainter. 
We first measure the magnitude difference Sm for each 
pair of overlapping images for each star in the matched 
astrometric catalogs. If we assume that this magnitude 
difference results entirely from systematic errors, then we 
can map out, in pixel coordinates, the magnitude correc- 
tions required to minimize all the Sm residuals. We do 
this by spatially binning the Sm values based on the pixel 
positions of the two measurements of each star. Weight- 
ing using the geometric mean of the errors in the two 
magnitude measurements, and rejecting > 3<r outliers, 
we find the maximum-likelihood magnitude corrections 
on a 100 x 100 grid in pixel space. 

The resulting maps for Bricks 15 and 21 are shown in 
Figure fl9] (top and bottom row, respectively) for F475W 
(left) ana F814W (right). If there were no spatially- 
varying systematics, the 5m values in each bin would 



be randomly distributed around zero and our correc- 
tion map would have small random values. Instead, we 
find coherent spatial patterns, at the ~ 0.03 mag level, 
that are consistent between bricks. The amplitude of 
the residual structure is clearly worse for F475W than 
for F814W. We believe that this difference most likely 
reflects the better large scale flat fields that are avail- 
able for F814W, which was calibrated against ground- 
based globular cluster observations (ISR ACS 2002-008 
by J. Mack et al.; ISR ACS 2006-06 by R. Gilliland et 
al.). Even in F814W, however, there are small < 0.01 
mag systematics associated with the chip gap; hints of 
these also have been reported in the analysis in ISR ACS 
2003-10 by R. van der Marel (see Figure 2). The maps 
in Figure [19] also indicate larger systematic errors where 
the boundaries of individual images lie in the distortion- 
corrected multi-drizzled reference images; stars that fall 
on these image edges are more likely to have unreliable 
photometry, and thus it is not clear whether the offsets 
at those positions reflects flat-fielding errors, or the dif- 
ficulty in doing reliable photometry near the edge of a 
chip. Even the largest offsets, however, are still less than 
5%. 

In future releases we expect to incorporate new large 



We see this effect for crowding errors as well (Section 4.7.2). 



scale flats based upon maps like those in Figure [19 
Unfortunately, we cannot generate equivalent maps for 
WFC3, due to the lack of overlapping data that covers 
the whole field of view for each chip. 

5. RESULTS 

5.1. Color Magnitude Diagrams 

The stellar populations revealed by the PHAT imaging 
are best understood through CMDs produced from the 
photometric catalogs. In this section we present CMDs 
for a range of radii within M31. To facilitate interpreta- 
tion of these data in the following discussion, we include 
in Figure [20] artificial CMDs generated for the PHAT fil- 
ter set, showing the positions of stars as function of their 
age and metallicity (for Ay = and a constant star 
formation rate, assuming photometric errors and limits 
based on artificial star tests from Brick 23). We also 
show the expected Galactic foreground populations in 
Figure [2TJ ca lculated with TRIL EGAL's default Milky 
Way model ( |Girardi et al. 2005[) (assuming no photo- 
metric errors, but applying depth cuts to approximate 
those in Figure 20). 

In what follows, we show CMDs for stars in Bricks 1, 
9, 15, and 21, drawn from the gst catalogs. For refer- 
ence, Bricks 9, 15, and 21 are sampling a series of M31's 
star forming rings at progressively larger radii, whereas 
Brick 1 is almost entirely dominated by the bulge (Fig- 
ure [I]) . The stellar surface density varies dramatically 
among these regions, with higher numbered bricks having 
systematically lower stellar surface densities, and thus 
deeper CMDs in the optical and NIR, due to reduced 
crowding (Figure 11). There are also significant varia- 
tions in the stellaraensity within Brick 1 itself, leading 
to a highly non-uniform limiting magnitudes within this 
single brick. 

5.1.1. Brick- Wide CMDs 

In Figures [22||25l we show Hess diagrams for the UV, 
optical, and NIRTfor stars in Bricks 1, 9, 15, and 21 (top 
left to bottom right). 
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The UV disk CMDs (Figure [22]) are dominated by a 
clear diagonal sequence of main sequence stars, extending 
down to F336W ~ 24.5 (corresponding to a ~3M star 
at the distance of M31). These young stars have ages 
of < 500 Myr (Figure 20 top right), assuming roughly 
solar metallicity for the young stellar populations. The 
main sequence starts to become more vertical brighter 
than F336W ~ 19.5 (> 15 M for solar metallicity main 
sequence stars) and has a color of F275W — F336W ~ 
—0.45 along the blue edge of the sequence. This color 
is slightly redder than the color of F275W - F336W ~ 



-0.6 seen in the unreddened simulations in Figure 20 



due to extinction from a Galactic foreground and dust 
associated with M31's star forming regions. 

The observed UV CMDs show a notable absence of 
stars redward of the main sequence, beyond the modest 
numbers expected from the Galactic foreground (Fig- 
ure |2l|). As can be seen in Figure 20 , even modestly 



sub-solar stellar populations can be expected to host a 
dramatic second sequence of blue core Helium burning 
stars, in addition to the main sequence. Additional sim- 
ulations (not shown) suggest that there should be a clear 
sequence of blue core Helium burning stars present for 
metallicities as high as [Fe/H]= —0.25. The absence of 
such a sequence suggests that even the outer disk of M31 
has close to solar metallicity at the present day. 

In contrast to the main sequence dominated CMDs for 
Bricks 9, 15, and 21, the bulge population in Brick 1 
shows a radically different UV CMD morphology. Al- 
though a smattering of main sequence stars are present, 
the most noticeable sequence is much bluer and more 
vertical, and dominates at colors of F275W — F336W ~ 
—0.5, fainter than F336W ~ 20. This sequence is made 
up of hot post-horizontal branch stars that had anoma- 
lously high mass loss rates on the RGB, which then al- 
tered what would otherwise be normal AGB evolution. 
The properti es of this popula tion are discussed exten- 
sively in |Rosenfield et al.| ( [2Q12| ) for the PHAT dataset. 



The optical disk CMDs in Figure [23] shows a much 
wider range of stellar ages than the UV (see Figure [20]) . 
Outside of the bulge (i.e., Bricks 9, 15, and 21), there 
is an obvious main sequence, extending nearly vertically 
at F475W - F814W - 0. The faintest detectable stars 
in the main sequence are found in Brick 21, and have 
typical masses of ~ 1.5 M , assuming solar metallicity. 
At redder colors, the stellar populations are dominated 
by a prominent red giant branch, punctuated by a well- 
populated red clump at F814W ^ 24.5 and a n otice- 
able AGB bump at F814W - 23.5 ( |Gallart|[l998| >. The 
RGB, red clump, and AGB bump are all smeared out 
diagonally by differential reddening. This effect is par- 
ticularly severe in Brick 15, which covers the main star- 
forming ring. The RGB terminates at bright magnitudes, 
with the tip of the red giant branch (TRGB) falling at 
F814W - 20.5 for bluer RGB stars, and F814W - 21.5 
for redder RGB stars. There is a significant population of 
AGB stars brightwards of the TRGB, with the bulk being 
no more than 1.5 magnitude brighter than the TRGB. 

There are also subtle, less populated features that can 
be discerned in the optical CMDs. First is a hint of a ver- 
tical sequence emerging brightward of the red clump and 
bluewards of the RGB. This feature is most likely due 
to core helium burning stars, which have characteristic 
ages of less than several hundred Myrs (see, for example, 



Fig. 1 of |McQuinn et al.|20lT ). These stars typically ap- 
pear as two narrow sequences in lower metallicity dwar f 
galaxies (see example CMDs in |Dalcanton et aL 2009), 
but in M31 the two sequences are redder and less dis- 
tinct, due both to higher metallicity and differential red- 
dening. As a result, the "sequence" manifests more as a 
subtle "edge" at F475W - F814W - 1.5 just above the 
red clump, corresponding to the blue edge of the core 
Helium burning evolution (see the solar metallicity pan- 
els of Figure 20 ) . The brighter stars in the core Helium 
burning sequences (i.e., classical supergiants) are difficult 
to separate from the Galactic foreground, which overlaps 
their positions in t he optical CMD (Figure [2l] see also 
Drout et al. (2009)), making it challenging to trace this 
feature to brighter magnitudes. 

In addition to the core Helium burning sequences, 
there is a more subtle secondary overdensity of stars ex- 
tending 0.5 mag faintward from the red clump, which 
is most evident in Brick 21 because of the lower level 
of crowding and better contrast in the CMD. This fea- 
ture corresponds to the "secondary clump" ( Girardi 



1999) caused by He-burning stars which are just massive 



enough to have ignited He in non-degenerate conditions 
at ages of 1 Gyr. 

The optical bulge CMD in Brick 1 shows some similar- 
ities with the disk CMDs. The bulge region is dominated 
by a very strong red giant branch with a clear TRGB, 
and a substantial population of AGB stars populating the 
region brighter than the TRGB. The RGB is also very 
broad, which is most likely due to a significant spread in 
stellar metallicity, given the lack of substantial differ en - 
tial reddening in the bulge (e.g., Montalto et al.||2009|) 



The crowding limit varies dramatically within different 
regions of Brick 1, however, such that the faint red clump 
is only detectable in the fields that are far from the center 
of the bulge. 

The bluer optical populations in Brick 1 differ sub- 
stantially from those seen in the disk fields, in two ma- 
jor ways. First, there is a noticeable diagonal sequence 
extending blueward from the RGB at F814W - 24, 
down to fainter magnitudes. In Figure [24] we overlay the 
post- AGB (P-AGB), post-early AGB (PE-AGB), AGB 
manque, and zero age horizontal branch evolutionary 
tracks that were used to analyze the UV populations in 
iRosenneld et al.| ( ^O^ p 7 ] The dashed portion of the 
cyan line indicates a veryfast interval of AGB manque 
evolution. Based on these tracks, it seems most likely 
that the diagonal sequence is made up primarily of blue 
horizontal branch stars (which are long-lived), though 
with a possible small contribution of AGB manque. How- 
ever, an exact attribution will depend on fuller model- 
ing that takes evolutionary lifetimes and likely reddening 
into account. We note that the stars in the diagonal se- 
quence are primarily found far from the center of the 
bulge, but it is not yet clear if this trend is due only to 
the increased depth at large radii, or reflects a real pop- 
ulation gradient; this populat ion is being analyzed more 
fully in (Williams et"at] fl2012| ) . 

The second major airier erice between the bulge and 
disk's blue populations is the behavior of the blue edge 
of the optical CMDs. In the disk regions, the upper blue 

27 PE-AGB and AGB manque stars are stars which lost too 
much mass on the RGB to undergo normal AGB evolution. 
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edge is defined by the nearly vertical main sequence. In 
Brick 1, however, there are few massive main sequence 
stars. Instead, the blue edge is curved to redder colors 
at bright magnitudes. Comparing to the tracks in Fig- 
ure [24j it appears that the blue edge in the bulge is de- 
finedby a combination of post- AGB and post-early AGB 
stars, with a possible contribution of main sequence stars 
at fainter magnitudes. 

Figure [25] shows the NIR color magnitude diagrams for 
all four bricks. The stars in these CMDs fall primarily 
on the older, redder RGB and AGB sequences, with an 
additional small contribution from the brightest main se- 
quence stars. The RGB is quite narrow in all the disk 
fields, and a clear TRGB is present at F160W - 18.3. 
However, differential reddening is still apparent, even at 
these longer wavelengths. The effect is most easily seen 
in Brick 15, where large amounts of dust lead the RGB 
to appear doubled, due to stars seen in front of and be- 
hind dust in the mid-plane of the galaxy. In the crowded 
bulge data, the RGB is broadened and does not extend 
more than ~3 magnitudes below the TRGB. In the outer 
disk, however, the CMDs are sufficiently deep that we 
see both a well-populated red clump and AGB bump. 
These features are extended to redder, fainter colors, 
suggesting that there is significant differential reddening, 
even at these large radii. Red core heli um burning stars 
are pr esent as well, but as discussed in |Dalcanton et aL 
(2012), they are difficult to distinguish from the main 
RGB and AGB sequences in this filter combination. 

The PHAT bricks are sufficiently large to contain a 
sizable number of foreground MW stars. These appear 
more numerous in the NIR CMDs, and appear as a nar- 
row vertical feature at F110W — F160W ~ 0.7 mag. Such 
a feature is evide nt in all our NIR C MDs. Simulations 
with TRILEGAL flGirardi et al.|2005[ ) reproduce well the 
shape and stellar density ol this feature (Figure 21). 



5.2. Global Structure of the MSI Disk 

Our resolved stellar data enable us to study the spa- 
tial distribution of selected stellar populations. Here we 
present initial results on the spatial distribution of the 
abundant metal-rich ([Fe/H] > —0.7) RGB stars with 
ages > 1 Gyr. Previous studies of the structure of M31's 
inner disk have relied on integrated light measurements 
of the combined stel lar population. Recent work by 
|Courteau et al.| ( |2011| ) uses / band and Spitzer 3.6 /im 
integrated light to fit the inner profile of M31. However, 
integrated light is an imperfect tracer of stellar mass, 
mixing stellar populations with varying mass-to-light ra- 
tios. While the NIR wavelengths used in recent work 
minimize the effects of dust, variations in M/L ratio at 
NIR wavelengths from young populations are still large 
due t o contribution from supergiant and TP- AGB stars 
(e.g., [Melbourne et al.|20"l2| . By using; star counts rather 
than integrated light we can reduce the contamination of 
these younger populations and directly study the struc- 
ture of an older, more dynamically mixed population. 

We select bright RGB stars using the NIR F110W and 
F160W filters, as this enables detection of even highly 
reddened sources. Figure [26] shows a NIR CMD of stars 
in Bricks 8 and 9, with our metal-rich RGB selection re- 
gion overlaid as a blue polygon. The blue edge of this 
region was selected to contain old RGB stars (~10 Gyr) 
up to ~1 magnitude below the RGB tip with metallicities 



[Fe/H] > —0.7. Younger RGB stars also fall within the 
box, with the solar metallicity 1.5 Gyr isochrone falling 
along the blue edge of the selection region. The redward 
extension of this box follows the reddening vector, and 
contains stars with E(B — V) < 3. Only the bright- 
est RGB stars were considered so that our selection box 
remains well above the completeness limit over all re- 
gions except the bulge. The RGB selection box excludes 
younger populations such as unreddened red supergiants 
and TP- AGB stars. 

After selecting the stars based on their position in the 
CMD, we calculated the densities of these stars in 30" x 
30" bins over all the regions with available WFC3/IR 
data. The RGB stellar densities are shown in the top 
panel of Figure 27 and range from >1 stars arcsecond -2 
in the bulge to ~ 0.01 stars csecond -2 in Brick 23. As 
would be expected from an old population, the RGB den- 
sity distribution is quite smooth, despite the outer fields 
covering several regions with strong star formation or 
high dust extinction. The axes in Figure [27] are rotated 
to follow the major and minor ax es assuming a position 
angle of 35° ( |Paturel et aT1|2003| ). 

We fit the 2-D RGB density of the stellar disk to an 
inclined exponential disk model after applying a small 
correction for dust extinction discussed below. We do 
not attempt to fit the bulge, due to incompleteness, and 
exclude the central region from our fits. Initial fits to the 
full data yield a disk scale length of 6.1 kpc. However, 
this fit is clearly biased by the presence of an overden- 
sity of RGB stars between 9 and 12.5 kpc, resulting in 
a model that significantly over-predicts the number of 
stars in Brick 21-23 (beyond 1 degree along the major 
axis). Excluding the region between 9.2 and 12.5 kpc, we 
obtain much better fits to the outer part of the galaxy, 
and a scale length of 5.3 kp c, in excellent ag r eemen t with 
the 5.3±0.5 kpc found by |Courteau et al.| ( [20TT] ). The 
best-fit position angle of the disk is 43.2°, with an incli- 
nation of 70°. This inclination is somewhat lower than 



the 74° inclination found by Courteau et al. ( 2011[ ) and 
the canonical value of 77° ass umed in earlier work (e.g., 
Walterbos fc Kennicutt| T988), perhaps due to our sam- 



pling an intrinsically thicker distribution of stars or our 
incomplete spatial covera ge. The posit i on an gle agrees 
well with that quoted by Sei gar et al.| ( 2008) 45°) and 
is slig htly rotated from the value given in |Tempel et al. 
( 2010[ 38.1°). The best-fit model and residuals are shown 
in the middle and bottom panels of Figure |27| The areas 
fitted are represented reasonably well by a single expo- 
nential disk: the reduced x 2 of this fit is 1.7, assuming 
Poisson statistics for the RGB star counts. In contrast, 
inclusion of the 10 kpc ring region results in a best-fit 
reduced x 2 of 3.7. 

Examination of the residual map in the bottom panel 
of Figure [27] shows a clear and significant overdensity be- 
tween 9 and 12.5 kpc along the major axis, with ~40% 
more stars than expected in our best-fit exponential 
model. This overdensity is associated wi th the "10 kpc 
ring" seen prominently i n UV emission (|Thilker et al. 
2005| and in the mid-IR QBarmby et al.|20U6p . This ring 



hosts ample recent star formation, and visual inspec- 
tion of the PHAT data shows the expected population 
of young main sequence stars in this region. However, 
since the stars selected for density mapping have ages 
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>1 Gyr, the overdensity observed here implies that the 
lOkpc ring, originall y identified through evidence for re- 
cent star formation (Habing et al. 1984 Gordon et al. 
2006), is overdense in more than just recently formed 
stars. 

A potential complication with this interpretation is 
that the higher extinction associated with this region 
could redden additional stars into our selection box. In 
Figure |26j the green box indicates areas of the CMD 
from which stars could be reddened into the selection 
box. This region includes primarily metal-poor RGB 
stars, AGB stars, and a small number of young super- 
giants. In the fitted region, the median number of stars 
in this box (referred to hereafter as C b s ) is 27% of the 
number of stars in our selection box (hereafter R Q bs)- 
The upper envelope in C b s /R bs is about 40%. Thus 
replicating the ~40% overdensity seen in the 10 kpc ring 
would require all of these stars to be reddened into the 
metal-rich RGB box. This scenario is implausible, es- 
pecially given that the RGB overdensity extends over a 
wide range of dust extinction. We can attempt to use 
the observed C b s /R bs to make a correction for the ef- 
fects of dust extinction on R b s . Assuming an intrinsic 
C/R ratio of 40%, we can calculate R%ntrinsic based on 
the observed C b s /R bs : 



R 



intrinsic 



R, 



obs ' 



1 + C bs/R< 



obs 



1.4 



(1) 



This correction does create a significantly smoother 
RGB density distribution, reducing the scatter in the fit 
residuals by ~20% both in the fitted regions and in the 
10 kpc ring. This reduction in residual scatter suggests 
that this correction does help account for the effect of 
dust extinction on our number counts and therefore we 
have included it in our fits discussed above and shown 
in Figure 27 However, the correction does not lessen 



the significance of the 10 kpc ring overdensity, at least in 
part because this region has a wide range of C b s /Robs- 
It therefore appears that the observed overdensity does 
in fact reflect a significant increase in the number metal- 
rich RGB stars at radii between 9 and 12.5 kpc along the 
major axis. 

In summary, we find disk scale lengths in agreement 



with previous estimates and find a strong overdensity 
(~40% over the expected density) of RGB stars affili- 
ated with the "10 kpc ring". These fits suggest that this 
structure is long-lived, and not simply an enhancement 
of recent star formation due to a spiral arm. A much 
more detailed analysis of the M31 structure as a func- 
tion of stellar population will be possible once the full 
data set is in hand. 



6. CONCLUSIONS 

We have presented the survey design, observations, and 
data products for the Panchromatic Hubble Andromeda 
Treasury. We have demonstrated the superb data quality 
offered by this new legacy of HST observations. Given 
that the optical and NIR data are crowding-limited over 
essentially all of M31's disk, the data at these wave- 
lengths cannot be surpassed in depth until higher res- 
olution facilities become available. Our hope is that the 
combination of wide wavelength coverage, photometric 
quality, and broad areal coverage will make this data set 
a rich scientific resource for decades to come. 
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TABLE 1 
Timing of Brick Observations 



Brick 


PID 


WFC3 East & ACS West 


WFC3 West & ACS East 


M31 Components Covered 






Start 


End 


Start 


End 




1 


12058 


2010-12-14 


2010-12-26 a 


2010-07-21 


2010-07-25 


Bulge, including nucleus 


2 


12073 






2011-06-30 


2011-08-03 


Outer bulge, NE star-forming arm 


3 


12109 










Outer bulge, inner disk, inner star-forming arms 


A 

4 


i oi n7 
iziu / 










NE star-forming arm 





1 9D7zL 






2011-07-13 


2011-08-03 


Inner-disk star formation 


u 


12105 










NE star-forming arm 


7 


12113 










Inner-disk star formation 


Q 
o 








2011-07-26 


2011-08-07 


IN J2j o Icli -1(Ji lillllg dilil, llllclloc old! 1L)1 IIldlKJIl, (JUlcI UloK 

Inner-disk star formation; significant overlap with brick 11 


Q 


1 9fl^7 


2010-12-30 


2011-01-01 


2010-07-12 


2010-07-16 


i n 


12111 










NE star-forming arm; outer disk 


11 


12115 










Smnntli intPT-firm Hi^lc* cHcni fir* Pint ovprlan witln triple Q 

UlllUU Ijll Ullid CX± 111 vllolv, Ol^lllllv^CXllL) UVC1 Idp VVllill VJL Iv^-tv *J 


12 


12071 






2011-07-01 


2011-07-13 


NE star-forming arm; outer disk 


13 


12114 










Smooth, mainly inactive inter-arm disk 


14 


12072 






2011-07-22 


2011-08-12 


NE star-forming arms; outer arm star formation 


15 


12056 


2011-01-10 


2011-01-25 


2010-08-07 


2010-08-15 


Richest sample of main-arm star formation 


16 


12106 






2011-08-29 


2011-09-03 b 


Main and outer-arm star formation 


17 


12059 


2010-12-29 


2011-01-09 


2011-07-26 


2011-08-06 


End of main star-forming ring 


18 


12108 






2011-08-21 


2011-08-28 


Outer-disk star formation 


19 


12110 










Inactive region in outer disk 


20 


12112 










Some outer-disk star formation; transition to halo 


21 


12055 


2011-01-10 


2011-01-15 


2010-07-26 


2010-07-28 


Outer-disk region of strong star formation 


22 


12076 






2011-08-14 


2011-08-17 


Outside of disk; transition to halo 


23 


12070 


2011-01-26 


2011-02-09 b 


2011-08-08 


2011-08-14 


Outside of disk; transition to halo 



a One orbit covering B01-F03/F06 for ACSF814WyUVIS were completed in 2011-01-03, after a guide star failure during the previous year's 
observing season. 

a One pointing has not yet been observed due to lack of guide stars (B16-F14/F17 for ACS/WFC3) or currently has invalid data due to failed 
guiding (B21-F16 for ACS); these observations are expected to be finished in the winter 2011/2012 observing season. 



TABLE 2 
Two-orbit Exposure Sequence 



ACS 


Exp 


Ax, Ay 


WFC3 


Exp 


Ax, Ay 




Filter 


(s) 


(pixels) 


Filter 


(s) 


(pixels) 


Notes 


F814W 


15 




F336W 


550 




F814W guard 


F814W 


350 




F275W 


350 










Large 37 Pixel Dither 


in AY to Cover WIS Chip Gap 


F814W 


700/800 


0.00, 0.00 


F336W 


700/800 






F814W 


455/550 


1.50, 1.50 


F275W 


575/660 




ACS Nyquist dither 


F475W 


10 




F160W 


399 


0.00, 0.00 


F475W guard 


F475W 


600/700 


0.00, 0.00 


F110W 


699/799 




ACS Nyquist dither 


F475W 


370 


1.01, -1.50 


F160W 


699 


1.52, 1.42 


ACS+WFC3 Nyquist dither 


F475W 


370 


0.48, -2.52 


F160W 


399 


1.94, 1.56 


ACS+WFC3 Nyquist dither 


F475W 


370/470 


-1.50, 1.00 


F160W 


399/499 


1.50, -0.01 


ACS+WFC3 Nyquist dither 



Note. — The F814W and WFC/UVIS exposures are obtained in the first orbit. The F475W and 
WFC/IR exposures are obtained in the second orbit. The dithers take place prior to the nominal 
exposure for which they are listed. Dithers are in pixels for the given camera. These are only specified 
for the sets of images for which dithers are designed to achieve the proper sub-pixel sampling needed 
to generate Nyquist-sampled images. The origin of the dithers is specified by the first image in the 
sequence. Several exposures have two values for the exposure time. The longer times are for images 
obtained with the "winter" sequence, which allows for longer orbital windows. 
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TABLE 3 
Approximate Corners of BRiCKS a 



Brick 


PID 


RAi 


Deci 


RA 2 


Dec2 


RA 3 


Dec 3 


RA 4 


Dec4 


1 


12058 


1 n 
1U 


Q7QCQ 

.o i yoy 


A 1 
41 , 


.ZOOoZ 


1U 


.OoO 1 1 


A 1 
41 


Q/1 


1 n 
1U 


.0 1 OlU 


A 1 
41 


9/1 QC7 


1 n 
1U, 


Q1 QQO 

.oioyz 


A 1 
41 


. lozyu 


2 


12073 


1 1 
1 1 


1 9090 

. izuzy 


A 1 
41 


1 8^/IQ 

. io04y 


i n 

1U 


.o ( ( yu 


A 1 
41 


.zooyy 


1 n 

1U 


.oioyy 


A 1 
41 


. 1D4DO 


1 1 
1 1 


.uoyzi 


A 1 
41 


.UooZ 1 


3 


12109 


1U 


.yozio 


/1 1 


.oozyo 


1U 


. i uoyu 


A 1 
41 


AAA 1 
.44410 


1 n 

1U 


p. a sno 


A 1 
41 


1A 1 7*3 
.o41 1 


1 n 

1U 


8Q1 1 8 


A 1 
41 


.ZDUDl 


4 


12107 


1 1 
1 1 


. ±yoo i 


/LI 
41 , 


, zozoy 


i n 

1U 


.youyz 


/L1 
41 


Q«/L Q/L 
.00404 


1 n 

1U 


oonon 
.ooyoU 


/LI 
41 


9^907 


1 1 

1 1 , 


, lozou 


/LI 
41 


1 Sfl/L/L 
. loU44 


5 


12074 


1 1 
1 1 


.U44oZ 


/1 1 
41 


A c;Q77 
,40O ( 1 


i n 

1U 


oni QQ 
.oUloo 


A 1 
41 


^^1 1 


1 n 

1U 


7/1 H9Q 
. 1 4uzy 


A 1 
41 


/I *397ft 
,4oZ ( D 


1 n 

1U 


noqci 

.yoooi 


A 1 
41 


^1 A A 
.00144 


6 


12105 


1 1 
1 1 


.Zo040 


/1 1 
41 


17A ni 

,o ( 4U1 


1 1 


n/i qc;q 


A 1 
41 


/I ^Sfi 


1 n 

1U 


Q0990 


A 1 
41 


.0O0O0 


1 1 
1 1 


99/1 qc 

,zz4yo 


A 1 
41 


971 89 
.Z ( loZ 


7 


12113 


1 1 
1 1 


. lo 1 I 4 


/LI 
41 . 


,0401U 


i n 

1U 


oQ/inQ 
.oy4uo 


/L1 
41 


«9««/L 
.0ZD04 


1 n 

1U 


.00Z 1 O 


/LI 
41 


,OZ4o4 


1 1 

1 1 , 


,U I DoU 


/LI 
41 


/L/L98/L 
.44Zo4 


8 


12075 


1 1 
1 1 


.O I o4o 


/1 1 
41 . 


/I CKA 1 A 
,4D4 1 4 


1 1 


. loooy 


/1 1 
41 


.o4o < y 


1 1 
1 1 


.U l 00 1 


A 1 
41 


/I /I /I <^9 
•444DZ 


1 1 
1 1 , 


Q1 P.7K 


A 1 
41 


.OOZOU 


9 


12057 


1 1 
1 1 


.zzyou 


A 1 
41 


.OoOoZ 


1U 


QOC07 

.yooo 1 


A 1 
41 


71 7flfi 
.111 UD 


1 n 

1U 


Q9/1 

.yz4oo 


A 1 
41 


fil /I SI 
.O14ol 


1 1 
1 1 


1 £81 ^ 


A 1 
41 


cqqi n 
•OoolU 


10 


12111 


1 1 
1 1 


.4 I 41U 


/LI 
41 . 


,oo4yo 


1 1 
1 1 


.zou 1 y 


/L1 
41 


.do ( iy 


1 1 
1 1 


. looyo 


/LI 
41 


COS 


1 1 

1 1 , 


/LI 91 ^ 

,41Z Id 


/LI 
41 


.40Zo I 


11 


12115 


11 


.28791 


41. 


,71003 


11 


.02353 


41 


.74143 


10 


.99982 


41 


.63066 


11. 


,26404 


41 


.59933 


12 


12071 


11 


.54277 


41 


.65341 


11 


.29918 


41 


.73581 


11 


.23717 


41 


.63373 


11. 


,48060 


41 


.55137 


13 


12114 


11 


.31717 


41 


.81950 


11 


.05228 


41 


.85054 


11 


.02883 


41 


.73975 


11 


,29355 


41 


.70886 


14 


12072 


11 


.57892 


41. 


,78777 


11 


.31430 


41 


.81941 


11 


.29041 


41 


.70867 


11. 


,55487 


41 


.67718 


15 


12056 


11 


.37644 


41 


,92498 


11 


.11115 


41 


.95616 


11 


.08756 


41 


.84538 


11 


.35268 


41 


.81435 


16 


12106 


11 


.63914 


41 


,89310 


11 


.37413 


41 


.92490 


11 


.35010 


41 


.81416 


11 


.61495 


41 


.78253 


17 


12059 


11 


.48855 


42, 


,02423 


11 


.22291 


42 


.05568 


11 


.19909 


41 


.94492 


11. 


,46456 


41 


.91363 


18 


12108 


11 


.75272 


41, 


,99292 


11 


.48736 


42 


.02497 


11 


.46310 


41 


.91427 


11. 


,72830 


41 


.88238 


19 


12110 


11 


.62005 


42 


,12072 


11 


.35409 


42 


.15246 


11 


.33001 


42 


.04174 


11 


,59581 


42 


.01015 


20 


12112 


11 


.88344 


42, 


.08995 


11 


.61776 


42 


.12230 


11 


.59324 


42 


.01163 


11 


,85876 


41 


.97944 


21 


12055 


11 


.68571 


42, 


,22514 


11 


.41935 


42 


.25704 


11 


.39512 


42 


.14633 


11. 


,66131 


42 


.11460 


22 


12076 


11 


.95062 


42 


.19503 


11 


.68455 


42 


.22754 


11 


.65987 


42 


.11690 


11 


,92579 


42 


.08456 


23 


12070 


11 


.86064 


42. 


.31678 


11 


.59401 


42 


.34909 


11 


.56944 


42 


.23843 


11 


,83591 


42 


.20629 



a All coordinates given in J2000. Corners are given for the approximate limits of the WFC3/IR coverage; 
WFC3/UVIS and ACS/WFC images will extend beyond these limits (see Figure [5}. 
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TABLE 4 
Dolphot Processing Parameters 



Description Parameter WFC3/UVIS ACS/WFC WFC3/IR 



Jr notometry Aperture oize 


DA 

KApGZ" 


o 
o 


o 
o 


o 


lllllci OKy XXdCIlUo 




15 


15 


g 


Outer Sky Radius 






oo 


9n 


^-statistic Aperture Size 


Rchi 


^ n 

o.u 


9 n 

z.u 


1 ^ 


PVintnmpfrAA Tvnp 
J. iiuiuiiicui v i y uc 


P^FPhnt 


1 


1 


1 


Fit Skv? 


r i u ua y 


3 


3 


3 


Spacing For Sky Measurement 


SkipSky 


2 


2 


2 


Qi rr m n (,1 inni n cr ThVit" ^llr^r 

kJlgllld VVllJJJJlllg, -L \JL kJxVy 




2.25 


2.25 


2.25 


Second Pass Finding Stars 


S e c ondP els s 


9^ 


9^ 


9^ 


>1 q rr n l n cr A lornrit nm 

kJCCXl ^lllllg, -Tilg^l ill 11111 


OtJCLX ^111 


I 
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Find/Make Aperture Corrections? 


ApCor 
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1 


Force Type 1/2? 


Forcel 
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1 


Align Images? 


Align 


2 


2 


2 


Allow Cross Terms in Alignment? 


Rotate 


1 


1 


1 


Centroid Box Size 


RCentroid 


1 


1 


1 


Search Step for Position Iterations 


PosStep 


0.1 


0.1 


0.1 


Maximum Single-Step in Position Iterations 


dPosMax 


2.5 


2.5 


2.5 


Minimum Separation for Two Stars for Cleaning 


RCombine 


1.415 


1.415 


1.415 


PSF Size 


RPSF 


10 


10 


10 


Minimum S/N for PSF Parameter Fits 


SigPSF 


3.0 


3.0 


3.0 


Make PSF Residual Image? 


PSFres 


1 


1 


1 


Coordinate Offset 


Psfoff 


0.0 


0.0 


0.0 


WFC3 PSF library 


UseWFC3 


1 




1 


ACS PSF library 


UseACS 




1 





TABLE 5 
PSF & Aperture Corrections 



Camera 


Filter 


Chip 


n psf 


PSF Correction 




Aperture Correction 


WFC3/UVIS 


F275W 


1 


102 


-0.0271±0.0047 


147 


0.0584±0.0176 


WFC3/UVIS 


F275W 


2 


93 


-0.0307±0.0044 


125 


0.0561±0.0167 


WFC3/UVIS 


F336W 


1 


115 


-0.0267±0.0049 


199 


0.0188±0.0154 


WFC3/UVIS 


F336W 


2 


111 


-0.0304±0.0048 


199 


0.0187±0.0142 


ACS/WFC 


F475W 


1 


204 


0.0051±0.0027 


198 


0.0048±0.0082 


ACS/WFC 


F475W 


2 


202 


0.0064±0.0021 


198 


0.0059±0.0083 


ACS/WFC 


F814W 


1 


231 


0.0091±0.0013 


200 


0.0159±0.0057 


ACS/WFC 


F814W 


2 


233 


0.0057±0.0010 


200 


0.0161±0.0053 


WFC3/IR 


F110W 




241 


-0.0059±0.0014 


199 


0.0573±0.0057 


WFC3/IR 


F160W 




241 


0.0062±0.0015 


200 


0.0784±0.0064 



Note. — n ps / and n ap are the median numbers of stars used to calculate the PSF and 
aperture corrections, respectively. PSF and aperture corrections are listed as the median ± 
the semi-interquartile range. All statistics have excluded data in Brick 1. The statistics for 
the short guard exposures in ACS have also been excluded. 
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Fig. 1. — Location of PHA T "bricks", superimposed on a Spitzer 3.6/i image jBarmbv et al.||2006| ), 24/x image dGordon et 111 2006), 
GALEX FUV image ( |Thilker et al.|2005| ), and Westerbork HI image (Brinks & Shane 19841, shown from top to bottom, respectively. The 
numbering scheme for the bricks is shown in the top panel; odd numbers are along the major axis, starting with Brick 1 on the bulge. 
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10 4 

MA) 

Fig. 2. — The relative transmission of the PHAT filter set (upper panel) as a f unction of wavelength . The lower panel shows a number 
of stellar spect ra for compariso n, including a Vega A-star spectrum |Bohlin 2007] blue), a G7III g iant (Castelli & Kurucz 2004 green), a 
M5III giant jFluks et al.|l994| red), and a carbon star with a 3200K atmosphere ( |Loidl et al.|2001| recT^ 




Fig. 3. — Map of the sub-pixel dither and exposure sequence for the two orbit visits. Arrows show the dithers relative to the previous 
exposure. The text at each vertex indicates the the filter and duration of the primary and parallel exposures taken at that position. 
Exposure times are given for the shorter schedlOO summer obs ervin g sequences and the longer sched60 winter sequences (i.e., "700/800s" 
means 700s in summer and 800s in winter; see text in Section |3.3| . For Orbit 1, only the sub-pixel dither between the third and fourth 
primary exposures is shown, since there is no dither between the first and second exposure, and the dither between the second and third 
is a large ~37 pixel dither to span the chip gap. 



Panchromatic Hubble Andromeda Treasury 



29 



Q 
>- 



0.6 



0.4 



0.2 



-0.2 





1 








-* 


. ~- 




~T~ ; ' 


_ F475W 

i i i 1 i i 


i 1 


i 1 i i i 1 i i 


-0.2 


0.2 


0.4 0.6 



X Dither (pixel) 

Fig. 4. — The sub-pixel dithers requested (red crosses) and achieved (points) are shown for ACS/WFC F475W images obtained during 
the Winter 2010 observing season. Note that the requested dithers do not fall exactly at the optimal 2x2 0.5 pixel locations, as the dither 
sequence operates on both the F475W and F160W niters simultaneously, and thus is a compromise between the needs of the two niters. 
The deviation of the actual dithers from the idea l interlace pa ttern is small enough that algorithms that can reconstruct Nyquist-sampled 
images from non-optimal dithers, such as that of |Lauer| ( |1999| >, should still work satisfactorily. 




Fig. 5. Maps of the relative exposure times across a brick for the WFC3/UVIS (left), ACS/WFC (middle), and WFC3/IR (right) 
tilings. The labeling of PHAT "fields" within a brick is shown superimposed on the WFC3/IR tiling. The left and right 3x3 half-bricks 
are taken with orientations 180° apart, as can be seen by the relative position of the dark circle in the WFC3/IR exposure map. The 
WFC3/UVIS exposures are dithered to cover the chip gap. The WFC chip gap in ACS is not dithered within a single orbit, but is covered 
by overlapping exposures. 
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Fig. 6. Four-color composite (F336W + F475W + F814W + F160W) images of Bricks 1, 9, 15, and 21 (upper left to lower right, 
respectively). The flux scaling is identical for Bricks 9, 15, and 21, to better show the radial change in the stellar density. Dark diagonal 
lines in the upper right and lower left corners are produced by the ACS/WFC chip gap (see Figure pi. Images were created by Zolt Levay 
(STScI). J ' 




Fig. 7. — Close-up of a false-color image of Brick 9, showing the rich detail available at full HST resolution. Courtesy Zolt Levay (STScI). 
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Fig. 8. Current PHAT coverage, as of Fall 2011, superimposed on a GALEX FUV+NUV image. Bricks 1, 9, 15, 17, 21, and 23 are 
largely complete, as are the first halves of Bricks 2, 5, 8, 12, 14, 16, 18, and 22. There are missing data in Brick 16 (Fields 14 and 17) and 
Brick 23 (Fields 13 and 16) that will be filled in during future observation cycles. 
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Fig. 9. — PSF corrections in magnitudes for all images in the survey to date, excluding the short guard exposures in F475W and F814W. 
For the WFC3/UVIS and ACS/WFC data, crosses and open circles indicate the corrections for chips 1 and 2, respectively. Median values 
are given in Table [5] The PSF corrections show no dependence on time or local stellar density. 
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Fig. 10. — Aperture corrections in magnitudes for all images in the survey to date, excluding the short guard exposures in F475W and 
F814W. For the WFC3/UVIS and ACS/WFC data, crosses and open circles indicate the corrections for chips 1 and 2, respectively. The 
aperture corrections vary with brick number, due to variations in the local stellar density, which decreases for increasing brick number, and 
for even-numbered bricks. The scatter tends to be systematically larger for bluer filters. The scatter is largest in Brick 1, which has the 
largest range of stellar densities internal to the brick, due to the presence of the center of the bulge in Field 10. Median values are given 
in Table [5] excluding both the guard exposures and images from Brick 1. 
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Fig. 11. — The magnitude of the 50% completeness limit, plotted as a function of radius for the high photometric quality gst catalog 
detections, for all fields observed with the WFC3/UVIS (horizontal blue points), ACS/WFC (lowest two purple sequences) and WFC3/IR 
(highest red sequences) cameras. The two filters in each camera are color-coded as in Figure [l2| The UV observations are photon-limited 
across the disk, with the exception of the very innermost regions of the central bulge field, and therefore have an essentially constant depth 
across the survey area. The optical observations are crowding-limited over most of the disk, and thus the depth of the ACS observations 
is a strong function of radius. The optical data does become photon-limited at a projected radius of ~0.9°, beyond which the depth of 
the observations becomes constant. The NIR observations are crowding-limited across the entire observed area, such that they are ~5.5 
magnitudes deeper in the uncrowded outer regions than in the bulge, in spite of having identical exposure times at all radii. 
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Fig. 12. — The magnitude of the 50% completeness limit as a function of the stellar surface density for the st catalog detections (S/N>4 
in a single filter; open circles), and the high photometric quality gst catalog detections (S/N>4 in both filters for a given camera, with 
additional stringent quality cuts; solid circles), for all individual fields observed with the WFC3/UVIS (left), ACS/WFC (center) and 
WFC3/IR (right) cameras by the end of the summer 2011 observing season (Table The two filters in each camera are color-coded as 
indicated in the individual panels. At low stellar surface densities, the limiting magnitude becomes constant, indicating that the depth is 
photon-limited and is independent of stellar density (seen in the UVIS and WFC plots). At high enough stellar surface densities, however, 
the number of detected (open circles) and/or well- measured (filled circles) stars becomes essentially constant, with the limiting magnitude 
being set by the point in the integrated luminosity function at which that maximum number of detectable stars is reached (seen in the 
IR and WFC plots); slight deviations from the constant stellar density are due to rapid changes in the stellar luminosity function at the 
TRGB and the red clump. The UV data (left panel) are clearly photon-limited (i.e., distributed horizontally on the plot), and have two 
different limiting magnitudes, because of the longer exposure times used in the winter observing season. The optical data (central panel) 
are photon-limited at stellar surface densities less than < 6 stars arcsecond - 2 , but are crowding-limited at most positions. The NIR data 
(right panel) never reach the theoretical photon-counting limit, and are crowding-limited everywhere. The central and right hand plots also 
show how the fraction of stars with well-measured photometry increases as the stellar density drops and the images become less crowded. 
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Fig. 13. — Median photometric errors reported by DOLPHOT in the six PHAT filters, for Bricks 1, 9, 15, and 21 (upper left to lower right). 
Medians are calculated for groups of 25, 400, or 200 stars, for the UV, optical, and NIR channels, respectively. Photometry is for stars 
with S/N>4 in a given filter, restricted to stars with well-measured photometry (FLAG=0) that do not fall near the chip gap or the edges of 
the image. Note that the formal photo metric uncertainties underestimate the true uncertainty, which is frequently dominated by crowding 
errors at the faint end (Figures |14)16| , and by systematic errors in the bright end; the horizontal dashed line indicates a characteristic 
amplitude of systematic uncertainties (Figure |18|). The steps at bright magnitudes are due to the limited number of significant figures used 
for the DOLPHOT photometric uncertainties. 
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Output Mag - Input Mag Flux Difference (Output - Input) / Flux Error Output Color - Input Color 

Fig. 14. — Comparisons between the input and output properties of artificial stars inserted into WFC3/UVIS images in Brick 1, Field 
9, for F275W (top row) and F336W (bottom row). Panels show the cumulative distribution of: (left) the difference between the recovered 
and inserted magnitude, such that negative numbers indicates stars that are recovered brighter than their true magnitude; (center) the 
fractional difference in flux, scaled by the photometric uncertainty reported by DOLPHOT such that positive numbers indicate stars that 
are recovered brighter than their true magnitude; and (right) the difference between the recovered and the input color of the artificial stars, 
such that positive values indicate stars that are recovered redder than their true color. Each curve shows the distribution for artificial stars 
with recovered magnitudes in a 1 magnitude wide range, spanning the total magnitude interval specified on each panel, color coded such that 
fainter stars are plotted with redder colors and thicker lines (i.e., the top row analyzes bins starting with stars having 20 < mp275W, out < 21 
(thin purple line), down to 24 < rap275W, o-ut < 25 (thick red line)). The center panel also includes a black line showing the expectation 
for a Gaussian with a width set equal to the median photometric error of stars in the bin. The UV data follow the Gaussian expectation 
well, because of the lack of crowding in the UV data everywhere outside the very center of the bulge. 
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Fig. 15. — Comparisons between the input and output properties of artificial stars inserted into ACS/WFC images in Field 9 of Bricks 
1, 9, 15, and 23 (top to bottom rows, respectively), for F475W (left column) and F814W (right column). Panels show the cumulative 
distribution of: (left) the difference between the recovered and inserted magnitude, such that negative numbers indicates stars that are 
recovered brighter than their true magnitude; (center) the fractional difference in flux, scaled by the photometric uncertainty reported by 
DOLPHOT such that positive numbers indicate stars that are recovered brighter than their true magnitude; and (right) the difference 
between the recovered and the input color of the artificial stars, such that positive values indi cate stars that are recovered redder than 
their true color. Line types indicate different bins in recovered magnitude, as described in Figure [14] Unlike the UV data, the optical data 
depart dramatically from a Gaussian due to the effects of crowding, in all but the outermost regions (see middle panel of each 3-panel plot). 
There is a tail of stars recovered with a flux that is significantly larger than their true flux, due to contributions from nearby stars. This 
bias can be several times larger than the DOLPHOT photometric error for faint stars in highly crowded regions. 
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Fig. 16. — Comparisons between the input and output properties of artificial stars inserted into WFC3/IR images in Field 9 of Bricks 
1, 9, 15, and 23 (top to bottom rows, respectively), for F110W (left column) and F160W (right column). Panels show the cumulative 
distribution of: (left) the difference between the recovered and inserted magnitude, such that negative numbers indicates stars that are 
recovered brighter than their true magnitude; (center) the fractional difference in flux, scaled by the photometric uncertainty reported by 
DOLPHOT such that positive numbers indicate stars that are recovered brighter than their true magnitude; and (right) the difference 
between the recovered and the input color of the artificial stars, such that positive values indi cate stars that are recovered redder than their 
true color. Line types indicate different bins in recovered magnitude, as described in Figure [l4] The WFC3/IR data depart dramatically 
from a Gaussian due to the effects of crowding, everywhere in the disk (see middle panel of each 3-panel plot). The effect is much worse 
than in the ACS data, due to the larger WFC3/IR pixel size. Moreover, the median star is recovered with a flux that is significantly fainter 
than its true flux, although crowding does induce the expected tail to brighter fluxes. 
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Fig. 17. — Top Row: Cumulative distributions of magnitude differences between repeated measurements of stars for Brick 15, binned into 
groups of stars with comparable magnitude (0.5 mag bin widths, with redder colors indicating groups of fainter stars), for WFC3/UVIS 
(left), ACS/WFC (center), WFC3/IR (right). Bottom row: lcr widths of the observed distributions of magnitude differences (black line, 
as inferred from the interquartile width to avoid skew from variable stars), as a function of the median stellar magnitude in a bin for Brick 
15. These widths are compared to the expected width from photometric uncertainties alone (blue line). The quadrature difference between 
them is an estimate of the systematic error, plotted as a red line. Systematic errors dominate the observed magnitude differences at bright 
magnitudes, and are larger at faint magnitudes, but still are typically small. The behavior of the magnitude differences is similar in other 
disk fields. Because the curves are based on pairs of overlapping stars, the plotted magnitude differences do not characterize the behavior 
of the full chip area for WFC3's UVIS and IR channels, which have limited spatial overlap. 




Fig. 18. — Inferred systematic errors as a function of magnitude for WFC3/UVIS (left), ACS/WFC (center), WFC3/IR (right), measured 
in Brick 15. The blue and red l ines indicate the systematic error in magnitude measurements for the bluer and redder filter in a single 
camera, respectively (see Figure [l7| . The black line shows the systematic error in the color. The systematic uncertainties are typically 
smaller in color than in magnitude, because many of the positionally dependent uncertainties (flat fields, PSF models) are likely to correlate 
between the two filters, and thus partially cancel out. The behavior of the magnitude differences is similar in other disk fields. Because the 
curves are based on pairs of overlapping stars, the plotted magnitude differences do not characterize the behavior of the full chip area for 
WFC3's UVIS and IR channels, which have limited spatial overlap (see Figure [5j. 
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Fig. 19. — Map of magnitude corrections as a function of position on the ACS reference image for F475W (left column) and F814W (right 
column), and Bricks 15 and 21 (top and bottom rows, respectively). The color bar on the right refers to the correction, in magnitudes, 
required to minimize the differences between our measurements for stars that appear in multiple images. The map of F814W offsets is 
quite smooth. In contrast, the F475W filter shows strong positionally-dependent variations that are consistent between the two bricks. 
These are likely to result from the lack of a large scale flat field correction for the F475W filter in the standard STScI flat field. 
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Fig. 20. Simulated CMDs for WFC3/UVIS (top row), ACS/WFC (middle row), and WFC3/IR (bottom row), for three different 
metallicities ([Fe/H] = —1.5, —0.75, and 0.0; left, middle, and right columns, respectively), with points color coded by age. The simulations 
assume a constant star formation rate, no foreground or internal extinction, and photometric errors and completeness derived from artificial 
star tests for outer disk fields in Brick 23. The solar metallicity CMDs are most likely to be characteristic of the young stellar populations 
throughout the disk, and the older stellar populations in the bulge. The inte rmediate metallicity ([Fe/H] = —0.75) CMDs are likely to be 
representative of the older thin and thick disk s tars (e.g., i Collins et al.|2011|), an d the low metallicity CMDs are likely to be characteristic 
of a potential stellar halo population (Chapman et al. 20U6 Kalirai et al. 2006a). Age color-coding is the same for all plots in single row, 
but varies from camera to camera. Vertical lines are plotted at a constant color for each c ame ra, and are shown for reference, to allow 
comparisons between different metallicities and with foreground populations shown in Figure [2l] 
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Fig. 21. Simulated CMDs for expected foreground populations in WFC3/UVIS (left), ACS/WFC (center), and WFC3/IR (right), 
calculated from TRILEGAL's default Milky Way model ( |Girardi et al.|[2005| >. CMDs assume a default area of 12 ; x 6.5 ; , corresponding 
to the 3x6 WFC3/IR footprint of a single brick; the number ot stars in each CMD is indicated in the plots. Stars in each panel are color 
coded by their distance modulus (m — M). The color and magnitude ranges are identical to those in Figure |20| Vertical lines are shown 
for reference, to allow comparisons with simulated stellar populations shown in Figure [20] 




Fig. 22. — UV Hess diagrams for complete Bricks 1, 9, 15, and 21 (from upper left to lower right). The Hess diagram is displayed with a 
square root scaling to highlight both high and low density features. Outside of Brick 1, the UV CMD morphology is dominate d by main 
seque nce stars. In the bulge-dominated Brick 1, however, a much bluer population due to hot HP-HB stars is also present (Rosenfield et al.| 
2012). The UV data reach a consist ent depth across all bricks, due to the low level of crowding. The vertical line is given tor reference, 
and matches those in Figures |20] and [2l| 
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Fig. 23. — Optical Hess diagrams for complete Bricks 1, 9, 15, and 21 (from upper left to lower right). The Hess diagram is displayed 
with a square root scaling to highlight both high and low density features. The optical CMDs show a varying mixture of red RGB stars and 
blue main sequence stars. The RGB also shows a clear red clump at F814W~24.5, and a brighter AGB bump at F814W~23.5. The depth 
in the optical depends strongly on radius because the ACS images are crowding-limited at most radii. The CMD for the bulge-dominated 
Brick 1 also shows a small diagonal sequence extending to F814W ~ 26, F475W — F814W ~ —0.5, due to blue horizontal branch stars 
( | Williams et al.|2012| >. The vertical line is given for reference, and matches those in Figures |20| and [2l] 
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Fig. 24. — Optical CMD of Brick 1, focusing on the diagonal sequence of blue stars. Overlayed are evolutionary tracks from Bressan et 
al. (in prep) for post-AGB, post early AGB, AGB manque, and the zero age horizontal branch stars (plotted as orange, green, light blue, 
and dark blue, respectively), assum ing zero foreground reddening with chemical composition Z = 0.07 and Y = 0.389 and an o:-enhanced 
composition typical of bulges (e.g., [Bensby et aL]|2010| ). We converted the tracks to the WFC3/UVIS photometric system following the 
|Girardi et al.||2008| > bolometric corrections updated with the latest WFC3/UVIS filter throughputs, and assumed m — M = 24.47. Stars 
are plotted as a Jless diagram in regions of high density on the CMD; howev er, to highlight the diagonal feature, the scaling of the Hess 
diagram and the transition to plotting points has been changed from in Figure [23] The dashed portion of the light blue AGB manque track 
indicates a very rapid phase of evolution; while this track appears to go through the overdensity of stars, the evolution through this phase 
is sufficiently rapid that the actual contribution of AGB manque stars to this feature is likely to be small. 
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Fig. 25. — NIR Hess diagrams for complete Bricks 1, 9, 15, and 21 (from upper left to lower right). The Hess diagram is displayed with 
a square root scaling to highlight both high and low density features. The NIR CMD is dominated by RGB stars, although AGB stars 
and MS are also present. The depth in the NIR depends strongly on radius because the NIR images are crowding-limited at all radii. In 
the outermost brick shown (Brick 21), the data reach below the NIR red clump (at F160W~23.5). There is also a clear AGB bump at 
F160W~22. Reddenin g an d extinction blur CMD features down and to the right. The vertical line is given for reference, and matches 
those in Figures [20] and [2l] 
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Fig. 26. — An example of selection of RGB stars from the NIR CMD. The blue selection region selects RGB stars with [Fe/H] > —0.7 
and reddenings up to E(B — V) = 3. This metal-rich RGB selection box only includes the brightest stars, to ensure uniform selection of 
RGB stars even in the crowded inner disk (i.e., Figure (ll). The green box indicates possible contaminants that could be reddened into 
the RGB selection area. These contaminants are metal-poor RGB stars, red supergiants and AGB stars. This box typically contains 25% 
of the stars in the metal-rich selection box. Red lines indicate RGB positions for [Fe/H] = (—2.3, —1.3, —0.7, 0.0) with an age of 10 Gyr. 
The NIR data are drawn from the merged catalog of Bricks 8 and 9; identical selection criteria were used in all other regions. A foreground 
extinction correction of E(B — V) = 0.062 was applied to the data before plotting. 
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Fig. 27. — Top: Surface density of luminous RGB stars falling in the cyan selection region from Figure [26] The overall RGB surface 
density distribution is quite smooth, with the exception of the inner region of the bulge, where incompleteness begins to affect even the 
brightest RGB stars. Middle: The best-fit single inclined exponential disk model, which was derived by fitting the observed surface density 
in the upper panel, after excluding the central regions of the bulge and the region form 9.2—12.5 kpc along the major axis. Bottom: Residual 
map showing the fractional deviation of the model from the data. Typical excursions are of order 10%, with the exception of the 10 kpc 
star forming ring, which appears as a true enhancement in the stellar surface density of ~40%, and is therefore a dynamical structure, 
rather than simply a localized enhancement in the rate of star formation. A literature value for the position angle of 35° was assumed in 
rotating the axes to be roughly along the major and minor axis. The fit along the major axis suggests a position angle of 43.2°. 



